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RESEARCH SUMMARY 


Background 

On October 31 » 1979» NASA launched a vector magnetometer into 
a fixed sun-earth orbit along the dawn-dusk meridian# at an 
altitude varying from 350 to 550 km. Although the mission's 
primary objective was to obtain data having sufficient accuracy 
and coverage to study static magnetization in the lithosphere# it 
has been recognized that there are often strong contributions 
from time varying external current sources in the magnetosphere, 
and perhaps the ionosphere# during magnetically disturbed 
periods. Hence this experiment presents an opportunity for 
studying external source fields as well as for studying the 
effects of electromagnetic induction in the solid earth. The 
purpose of this report is to describe progress in employing 
satellite observations to characterize the internal contribution 
to these rapidly varying# transient field effects - the magnetic 
field due to time-v.irying electromagnetically induced currents 
flowing in the earth. 

Theoretical Assessment Induction Effects Al Satellite 
Altitudes 

Hermance (1982) considered from a theoretical viewpoint 
several questions related to satellite-induction studies. First . 
are# in fact, the magnetic field effects of induced currents 
flowing in the earth significant for near-earth satellite 
observations? He found that for distant sources in the 



magnetosphere r magnetic fields from electric currents induced in 
a spherical earth may contribute from 30 to 40% of the external 
source field amplitude. For near-earth sources beneath the 
satellite^ e.g. in the ionospherer the external and internal 
fields tend to cancel r resulting in a much smaller total 
horizontal field than would be associated with sources above the 
satellite^ or for the sources alone» if the contribution from 
induced currents in the earth were negligible. Second . Hermance 
(1982) investigated the effects at satellite altitudes of lateral 
differences in the gross conductivity structure of the earth. He 
found, for the particular case when the current sources are in 
the ionosphere (i.e., between the satellite and the ground), that 
local concentrations of induced current may generate unique 
magnetic field signatures at satellite altitudes. Although the 
absolute amplitude of the field components may be an order of 
magnitude smaller than for the case when the source is above the 
satellite, the magnetic field perturbation measured by the 
satellite as it crosses the geologic contact is significantly 
larger in relative terms for ionospheric sources. On the other 
hand, for the ca''' of remote magnetospher ic sources (e.g. the 
ring current fields) , he pointed out that the contribution from 
lateral heterogeneities are, for all practical purposes, 
negl igible . 

In short, the theoretical studies of Hermance (1982) showed 
that a) the effect? of induced fields were significant at near- 
earth satellite altitudes, b) the position of the source current 
(i.e. whether it was in the ionosphere, between the earth and the 
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satellitC/ or whether it was remote in the magnetosphere) had a 
profound effect on the electromagnetic response of the structure, 
and c) if, in fact, the source was remote (e.g. due to the 
magnetospheric ring current) , then the effects of lateral 
heterogeneities are essentially negligible. 

Qbseivatipnal £yl dance inx induction Effects in magsat onln 

Hermance and Rossen (1982) pursued these questions from an 
observational viewpoint. They identified a number of issues that 
needed to be addressed before the efficacy of satellite induction 
studies could be established. One was whether satellite data 
offered an improvement over ground-based observations. They 
found that the "noise" level during magnetic storms was 
significantly less at satellite altitudes, particularly on the Z 
component. Another issue they addressed was the presence of 
ionospheric current systems during major magnetic storms. In 
comparing ground-based and satellite observations of the Z and H 
components during a major magnetic storm during the MAGSAT epoch, 
they found little evidence for zonal ionospheric current systems. 
Although they detected a significant offset in the horizontal 
field when satellite and ground-based data were compared, which 
they originally felt indicated a significant contribution from 
ionospheric currents, they later traced this offset to an 
artifact introduced by uncertain base-lines. Present feeling is 
that there is little evidence that significant zonal ionospheric 
currents flow in the dawn and dusk sectors during the main and 
recovery phase of major magnetic storms. 
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Hermance and Rossen (1982) proceeded to calculate "primitive" 
electromagnetic response functions and determined an average 
radius of 5,370 (± 120) km for an equivalent super-conductor, 
representing the conducting region in the earth's interior. In 
addition they found little differences between the 
electromagnetic response of oceanic areas and other areas of the 
globe. This prompted them to more closely scrutinize MAGSAT data 
for the effects of lateral conductivity heterogeneities in the 
earth. 

To evaluate the effects of lateral heterogeneities, Hermance 
and Rossen (1982) investigated the well-known ocean-edge effect 
in the vicinity of Australia. They found that the effect is of 
such a restricted scale-size laterally, that it is damped to 
negligible amplitudes at the altitude of MAGSAT - this 
corroborates the conclusions that Hermance (1982) reached from 


theoretical considerations. 
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Conclusions 

The work summarized here is best described as the 
reconnaissance phase of evaluating the contribution that 
satellite observations can make to studying electromagnetic 
induction in the solid earth. Although certain elements of our 
analysis are qualitative# several points can nevertheless be 
made. First . satellite data apparently suffer far less from the 
effects of near surface lateral heterogeneities in the earth than 
do ground-based data. Second , ionospheric currents during the 
recovery phase of major magnetic storms appear to be miniral# at 
least in the dawn and dusk sectors where MAGSAT was flown. Hence 
the internal contributions that satellites observe during these 
times is in fact due primarily to induction in the earth with 
little or no contribution from ionospheric currents. Third . the 
interpretation of satellite data in terms of primitive 
electromagnetic response functions suggests that more 
quantitative interpretations hold great promise for refining 
models of the electrical conductivity in the upper mantle. Even 
the grossly over-simplified interpretation presented here results 
in a surprisingly well-resolved radius for an equivalent super- 
conductor representing the conductivity region of the eorth's 
interior (5,370 ± 120 km). This is indeed encouraging for 
follow-on work in this area. What is required now is a refined 
spectral analysis of a number of different transient events, 
accounting for a more complicaced source field than has been 
assumed hitherto. Having resolved the frequency dependence of 
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the electromagnetic response parameter » one will be able to 
investigate the global conductivity structure of the earth's 
interior using more realistic models than the one used here. 
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Summaiy of Attached Documentation 

Attached as Appendices to this Final Project Report are the 
following items; 

A ppendix 1: Abstracts of Oral Presentations at National and 

International Meetings. This is a collection of four 
oral presentations describing research of which a 
significant portion was supported by NASA funds. 

A ppendix 2.: Model Simulations of Possible Electromagnetic 

Induction Effects at MAGSAT Altitudes. A reprint of a 
paper which appeared in Geophvs . Res . i, 

373-376, 1982. 

A ppendix 1; Global and Regional Electromagnetic Induction 

Effects in MAGSAT Satellite Data. A revised version of 
a paper read at the 6th Workshop on Electromagnetic 
Induction, held in Victoria, B.C., during August of 
1982. 

Appendix a: Electromagnetic Induction by Finite Wave-Number 

Source Fields in 2-D Lateral Heterogeneities; The 
Transverse Electric Mode. A preprint of a paper 
submitted to Geophys . J.fi. 1983. The 

development of this theoretical effort was prompted by 
research supported by several funding agencies. 


REFERENCES 


Banks » R.J., Geomagnetic variations and the electrical 

conductivity of the upper mantle # G ^ophys . fioy. Aslion. 

Soc. . n, 457-487, 1969. 

Chapman, S. and J. Bartels, Geomagnetism . Clarendon Press, 

Oxford, 1940. 

Didwall, E.M., The electrical conductivity of the Earth's upper 
mantle as estimated from satellite measured magnetic field 
variations, Ph.D. thesis, Johns Hopkins University, Baltimore, 
MD, 1981. 

Didwall, E.M. The electrical conductivity of the upper mantle as 
estimated from satellite magnetic field data with implications 
for temperature. Submitted to Geophys. Res . . 1983. 

Hermance, J.F., Electromagnetic induction and source field 
effects in MAGSAT data (Abstract), Programme and AbStlfl£t£f 
Fourth lAGA Scientific Assembly, edited by N. Fukushima, p. 
538, lUGG Publ. Office, Paris, France, 1981a. 

Hermance, J.F., Are there induction effects in Dst corrections to 
MAGSAT data? (Abstract), 12, p. 269, 1981b. 

Hermance, J.F., Finite source fields coupled to lateral 

conductivity heterogeneities; Effects on magnetotelluric and 
magnetic gradiometric deep sounding experiments (Abstract) , 
Programme and Abstracts . Fourth lASA Scientific Assembly r 
edited by N. Fukushima, p. 566, lUGG Publ. Office, Paris, 
France, 1981c. 

Hermance, J.F., Model simulations of possible electromagnetic 

induction effects at MAGSAT altitudes, Geophy s . Ban. Lnti., i, 
373-376, 1982. 

Hermance, J. F. and M. Rossen, Regional and global 

electromagnetic induction effects in MAGSAT satellite data 
(Abstract) , Sixth Workshop on Electromagnetic Induction in the 
Earth and Moon, Working Group I -S, lAGA, University of 
Victoria, British Columbia, Canada, 1982. 

Hermance, J.F., Electromagnetic induction by finite wave-number 
source fields in 2-D lateral heterogeneities; The transverse 
electric mode, Manuscript in preparation, 1983. 

Langel, R.A., Internal and external storm-time magnetic fields 
from spacecraft data, lAGA Program Ab.St.I,.fl,.C.tS iDX XYl Geneifll 
AaS-iilQtUj. , Genoble, France, 1975. 

Langel, R.A., R.H. Estes, G.D. Mead, E.B. Fabiano, and E.F. 

Lancaster, Initial geomaonetic field model from MAGSAT vector 
data, Res . Lett ■ . 2, 793-796, 19G0. 


Langelr R.r J. Herbert, T. Jennings, and R. Horner, MAGSAT 

processing; & report for Investigators , Technical Kcmorandum 
82610, 178 pp., NASA, 1981. 

Langel, R.A., The magnetic earth as seen from MAGSAT, initial 
results, Geophys . £££. LsH., 2., 239-242, 1982. 

Langel, R.A. , G. Ousley, J. Herbert, J. Murphy and M. Settle, The 
MAGSAT mission, Geophys . Res . Lett . . i, 243-245, 1582. 

Langel, R.A., J.D. Phillips, and R.J. Horner, Initial scalar 
magnetic anomaly map from MAGSAT, Geophys . Res . Lett . . 2., 
269-272, 19B2. 

Langel, R.A., C.C. Schnetzler, J.D. Phillips, and R.J. Horner, 
Initial vector magnetic anomaly map from MAGSAT, Geophys . Res . 
Lfilt, 2, 273-298, 1982. 



APPENDIX I 


Abstracts of Oral Presentations at 
National and International Meetings 




Abstract of paper presented at the Spring Meeting of the American Geophysical 
Union, 1981. 


ARE THERE INDUCTION EFFECTS IN Dst CORRECTIONS 
TO MAGSAT DATA? 

John F. Hermance (Department of Geological 
Sciences, Brown University, Providence, 

Rhode Island 02912) 

There is a need, when deriving residual 
field models from MAGSAT data, to be able to 
apply corrections to the data for the con- 
tribution from Dst. Unfortunately induction 
effects in Dst preclude in general the simple 
correction of the observed field by subtracting 
the Dst-term; cns needs to account for the in- 
duced phase-shifted internal component. In 
addition, the temporal behavior of Dst is such 
that a magnetic storm with a sudden commence- 
ment has higher frequency energy at early 
storm-time (T < 2 days) and lower frequency 
energy at late storm-time (2<T<5 days). The 
result is that the ratio of internal to exter- 
nal fi.elds is not a simple constant, but 'e 
plitude of the ratio changes with time during 
magnetic disturbance. Fortunately, simple 
model calculations indicate that during late 
storm-time the internal : external field ratio 
is essentially constant so that there is hope 
that quantitative correction factors can be 
applied to late storm-times in a straight- 
£(jrward way. As one attempts such a correction 
during earlier storm-time, the morphology of 
the Dst event has to be taken into account. 



Abstract of paper presented at Fourth Scientific 
Assembly of lAGA, Edinburgh, 1981. 

ELECTROMAGNETIC INDUCTION AND SOURCE FIELD EFFECTS IN MAGSAT DATA 


J.F. Hermance (Department of Geological Sciences, Brown 
University, Providence, RI 02912 


It is advantageous when compiling regional magnetic charts to be able 
to use satellite data during periods when the magnetic field is modestly to 
severely disturbed. It is necessary, therefore, to develop algorithms to 
compensate observed field data for the effects of time varying magnetic sources 
of external origin, particularly accounting for possible induction effects in 
the earth. 

To illustrate some of the complications in such a procedure, a series 
of theoretical models is presented to Illustrate the effects of external 
. jrce fields coupling to an earth having a finite conductivity. Among the 
cases considered are: 

1) Induction effects in the Dst index such that a storm-time dependent 
correction factor need be applied. 

2) The electromagnetic coupling of ionospheric current systems to both 
a homogeneous earth having finite conductivity, and to an earth 
having gross lateral variations in its conductivity structure, e.g. 
the ocean-land interface. 


Abstract of paper presented at Fourth Scientific 
Assembly of lAGA, Edinburgh, 1981. 


FINITE SOURCE FIELDS COUPLED TO LATERAL CONDUCTIVITY 
HETEROGENEITIES: EFFECTS ON MAGNETOTELLURIC AND MAGNETIC 

GRADIOMETRIC DEEP SOUNDING EXPERIMENTS 

J.F. Herman ce (Department of Geological Sciences, Brown 
University, Providence, Rhode Island 02912) 


Studies of the dynamical evolution of source fields in the ionosphere 
and magnetosphere are in many cases incompatible with global scale sources. 

In some cases, source fields have quadrant size dimensions, in others coherent 
scale distances of a few .mndred km, even at mid-latitudes. The possible 
effect of these finite source dimensions on induction experiments therefore 
require renewed scrutiny, particularly as deep-sounding measurements are 
performed at periods greater than 10^ sec. The magnetic variation gradio- 
metric technique appears to suffer from source et'’'?cts at least as much, 
and perhaps greater than, the magnetotelluric te--'-. ue. When finite source 
fields are coupled to lateral heterogeneities, the .iteraction becomes quite 
complex. We will report on our progress in generals iing the interaction 
between the spat J 'I wave-length of the source field and the scale-size of 
lateral heterogeneities. This is used to evaluate the conditions under 
which various phenomena are more dominant than others. 


Abstract of paper presented at the 6th Workshop on Electromagnetic 
Induction in the Earth and Moor, Victoria, British Columbia, 1982. 


Regional and Global Electromagnetic Induction Effects 
in MAGSAT Satellite Data 


John F. Hermance (Department of Geological Sciences, 
Brovn University, Providence, Rhode Island 02912) 

Michael Rossen (Department of Geological Sciences, 
Brown University, Providence, Rhode Island 02912) 

Our present analysis of MAGSAT magnetic data is 
separated into two stages. Both stages are designed 
to study the extent to which we can detect and con- 
strain the nature of electromagnetic induction ef- 
fects using satellite observations. 

The first stage consists of comparison of data 
from contiguous orbit passes over large scale geo- 
logic boundaries, such as ocean-land Interfaces, 
at several levels of magnetic disturbance. The 
purpose of these comparisons is to separate induc- 
tion effects from effects of lithospheric magnetiza- 
tion. The procedure for reducing the data in order 
to make such comparisons is the following: 1) iden- 

tifying and subtracting quiet time effects; 2) mod- 
elling and subtracting first order ring current 
effects; 3) projecting an orbit track onto a map 
as a nearly straight line so it can serve as an 
axis on which to plot the corresponding orbit pass 
data in the context of geography. 

The second stage of our analysis consists of 
comparison of MAGSAT data with standard observatory 
hourly data. The purpose is to constrain the time 
evolution of ionospheric and magnetospheric current 
systems. Qualitative features of our ground based 
dataset are discussed. The methods for reducing 
the ground based data consists of: 1) a baseline 

is removed from the data; 2) a spherical harmonic 
analysis is performed on the data; 3) the spherical 
harmonic coefficients are used to reconstruct, using 
a suitable method of interpolation, an equivalent 
"ground pass" dataset corresponding to each orbit 
pass dataset under study; 4) a separation of the 
reconstructed "ground pass" data into internal 
and external components is made; the data are then 
projected up to MAGSAT altitude with the assumption 
that there are no Intervening current systems. 

Comparisons between ground passes and orbit passes 
are used to differentiate between the presence of 
ionospheric and magnetospheric current systems. 

Special attention is payed to analysis of long 
period effects, such as the storm recovery phase, 
with which the inclusion of hourly values in our 
dataset is most effective. 
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GLOBAL AND REGIONAL ELECTROMAGNETIC INDUCTION EFFECTS 
IN MAGSAT SATELLITE DATA 


Johr F. Hermance 
Michael Rossen 

Geophysical, /Electromagnetics Laboratory 
Department of Geological Sciences 
Brown University 
Providence, Rhode Island 02912 


Revised version of a paper read at the 
Sixth Workshop on Electromagnetic Induction 
Victoria, British Columbia 
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INTRODUCTION 


On October 31# 1979# NASA launched a vector magnetometer into 
a fixed sun-earth orbit along the dawn-dusk meridian# at 
altitudes varyng between 350 km and 550 km. The resolution of 
the data is to within ±1 nT. The mission's primary objective was 
to obtain data having sufficient accuracy and coverage to produce 
a detailed map of static magnetization in the lithosphere (Langel 
et al. # 1982) . 

It has been recognized that time varying contributions from 
external current sources are present in the data. This presents 
an opportunity for earth induction studies using vector magnetic 
variation data - Langel (1975)# Didwall (1981# 1983)# and others 
have attempted such analyses ptex'iously using total field scalar 
data from earlier satellites. 

A number of issues need to be addressed before the efficacy of 
satellite induction studies can be established. First# do 
satellite measurements result in a data base which has a 
significantly lower noise component than ground-based 
measurements? Second# we must determine whether ionospheric 
currents are present at nominal altitudes of 110 km during 
magnetically disturbed periods. If present and unconstrained, 
such currents, since they flow beneath the satellite's orbit, 
would tend to bias satellite measurements of terrestrial fields. 


Finally, we must ask whether the effects of electromagnetic 
induction in lateral heterogeneities in the earth are significant 
at satellite altitudes. In this report, we describe our progress 
in resolving these issues. 

MAGNETIC STORM-TIME VARIATIONS 
Reduction q 1 Ground-Based Data 

For the purposes of our present analysis we want to compare as 
directly as possible ground-based data from standard magnetic 
observatories to satellite data. Ground data is readily 
available from the greatest number of observatories in the form 
of hourly averages. Although these represent time series which 
are somewhat courser than desirable, they seem to be satisfactory 
for the type of analysis presented here. Our ground data set, 
therefore, consists of standard observatory hourly means f jr 32 
sites between 55° and -55° geomagnetic latitude, as shown in 
Fig. 1. 

The data being examined are from the period between 
0:00 February 12 UT and 0;no February 22, 1980. This epoch 
includes one of the most well defined geomagnetic storms to occur 
during the MAGSAT mission, with sudden commencement occurring at 
3:09, February 14 UT (preliminary indent if ication ; Langel 
1981). In Figure 2 we see the H, D, and Z data from the Tucson 
standard magnetic observatory in the southwest United States. 

The H-component in particular shows a well-developed main phase 
on February 16-17, followed by a slow recovery phase to normal 
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ba8e**line values after a period of 5-6 days. Superimposed on the 
storm-time magnetic disturbance one can readily discern the 
periodic diurnal variation (Sq ) , which is of course strongest 
during local daylight hours. 


Our method of reducing the ground based data is as follows: 

1) All data are cast in terms of geomagnetic coordinates and 
geomagnetic field components (H, Df Z) . 

2) A baseline f calculated as an average of two hours about 
the quiet time local midnight immediately before the 
onset of the magnetic disturbance (see for example Figure 
2 ), is removed from each component at each observatory. 

3) The vector components at each observatory are 

geometrically corrected for the various latitudes at 
which they were recorded. For the lowest order surface 
harmonic, we assume H « H cos and Z « Z sin 

The coefficients H and Z are determined by dividing all 
H-fields by the cosine of the latitude a. id all Z-fields 
by the sine of the latitude. 

4) Time series are then generated for each component as a 

function of storm-time for local tine sectors: dusk, 

midnight, and dawn. This is done by calculating, for the 
latitudinally compensated data from each observatory, 
two-hour means centered at local times of 1800 hr (dusk), 
0000 hr (midnight) and 0600 hr (dawn) . Each observatory 
therefore contributes one sample in each of the three 
local time sectors for each day of storm-time. 


5) The time series for the dusk and dawn sectors are then 



directly compared with the reduced satellite data 
described below. 

The ground-based time series for the H-component recorded 
during the dusk, midnight and dawn sectors, respectively, are 
shown in Figure 3. The midnight samples are clearly the best 
behaved, but all three samples display well-developed classical 
m^in and recovery phases. MAGSAT data are limited by the nature 
of the orbit to dawn and dusk sectors, so we are constrained to 
using ground-based data at only these two local times. Of these 
two times, dawn samples seem to display a more consistent pattern 
during pre-storm times as well as during the recovery phase. The 
veil-known dawn-dusk asymmetry in main phase amplitudes is clear 
in Figure 3. 

Both the Z-amplitudes and the H-amplitudes for ground-based 
data during the dawn sector are shown in Figure 4. During pre- 
storm time the baseline for the vertical field is well-described, 
however during the recovery phase a greater scatter of points is 
apparent - we ascribe this to the effects of lateral 
neterogeneities within the earth which tend to affect the 
vertical field more strongly than the horizontal field. 

Reduction £Lf MAGSAT 

Five second averaged data, with the standard reference field 
model removed, were taken from Investigator-B tapes to form our 
initial dataset. Spikes were automatically culled from the data. 
A rotation into the geocentric magnetic dipole reference frame 
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was then applied to both the coordinates of the satellite 
position and to the components of the vector magnetic field. 
Following convention^ H« D and Z refer to the magnetic Northf 
magnetic East and vertical down field components respectively. 
These basic data were then reduced in the following way: 

1) H data were fit to a cosine function of magnetic latitude 
and Z data were fit to a sine function of magnetic 
latitude. The ascending (dusk) and descending (dawn) 
part of each orbit pass between -45 deg and 45 deg 
magnetic latitude were fit separately to allow 
observation of the dawn-dusk asymmetry. 

2) The amplitudes of these cosine and sine fits of 
succeeding passes were gathered together as new, smoothed 
datasets, directly comparable to the ground based field 
coefficients described above. The equator-crossing times 
of the passes were used to form a time base for the 
datasets. This provided a sampling interval of about 
46.5 min, or one half the orbit period. Two datasets 
were thus created - a dawn dataset for H and Z and a dusk 
dataset for H and Z. The geographic location of the 
actual passes used for this analysis are shown in 
Figure 5. 

At this stage a direct comparison of the reduced satellite and 
ground-based data sets (Figure 6 and 7) shows that while the 
general form of the appropriate components is quite similar, a 
persistent offset between the two datasets seems to be present. 
Using hourly Dst values from Langel £1 al. (1981) , a mean Dst 


value was computed for these two time intervals to determine 
whether an offset due to residual Dst field effects should be 
expected. The arithmetic mean of Dst indices for the pre-storm 
period (0:00 February 12 UT to 0:00 February 14 UT) was 
0.9 gammas. The mean for the two and a half days at the end of 
the recovery phase (12:00 February 19 UT to 0:00 February 22 UT) 
was -10.8 gammas. 

We conclude^ therefore, that the Dst index, particularly at 
the beginning of the storm is close to zero. Langel and Estes 
(1983), however, have found in an independent analysis that 
significant external field is present even when Dst is zero. 
Employing a potential field separation of the MAGSAT satellite 
data and performing a regression of the external potential term 
on the Dst index, they find a correlation of the form: 

= 19.0 - 0.65 Dst (nT) . 

It does not seem unreasonable, therefore, that offsets of 
10-20 nT may be present in both the ground-based and satellite 
data, due to long term quasi-state external field contributions. 
Not wishing to address the issue of establishing absolute 
baselines to either the ground data or the satellite data at this 
time, we can nevertheless calculate relativ.e offsets between the 
two data sets in the following way. 

Feeling that the satellite data is least "noisy" - except for 
the constant offset - a smooth curve (a cubic spline) was fit in 
a least m.ean square sense to the first two days of the pre-storm 
satellite data (the same times indicated above) . A similar curve 
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was fit to the last two and a half days of the recovery phase. 

The residual could then be calculated between actual ground-based 
data observed at arbitrary times during these two periods and 
satellite observations interpolated from the smooth cubic spline 
functions. 

The residuals for each component were sorted into bins 
corresponding to the prestorm period, the late recovery phase 
period, and the combination of both periolc. In addition, the 
data were sorted into dawn passes and dusk passes. Histograms of 
the residuals are reproduced in Figures 8-11. 

The average residuals from this analysis are summarized in 
Table 1. Their magnitudes and signs generally agree with the 
analysis by Langel and Estes mentioned above. The next and final 
step in reducing the satellite data therefore was the following. 

3) The offsets indicated in Table 1 were subtracted from the 
corresponding satellite components. An example of one of 
the resulting data sets for Z and H is presented in 
Figure 12. 

This represents the final reduced satellite data. Note that 
the daca from both components show significantly lower scatter 
than the corresponding ground-based data in Figure 4. This 
comparison graphically depicts the potential advantages in using 
satellites to study magnetic storm phenomena - particularly earth 
induction effects. We feel that a substantial contribution to 
the noise in the ground-based data of Figure 4 is due to effects 
of lateral conductivity heterogeneities within the earth. The 
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field romponents at some observatories may suffer from 
contamination due to edge-effects encountered near 
discontinuities in the electrical structure? in other cases, the 
fields may be perturbed by different conductivity structures 
beneath the observatories. 

The satellite data and ground-based data for the dawn and the 
dusk sectors, respectively, are actually overlaid in Figures 13 
and 14. These figures again illustrate the lower scatter in the 
reduced satellite components. However, in spite of the somewhat 
higher noise level in the ground data, there are specific time 
intervals and events for which the tracking of the two datasets 
is quite remarkable. In spite of the simplicity of our approach, 
the results are encouraging. 
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£2n ths. Possibility i.ongfipiifiii,c Currents 

We now address the issue as to whether appreciable ionospheric 
currents flow at low and mid-latitudes during disturbed periods. 
Ionospheric currents flow at a nominal altitude of 110 km aid 
during magnetically quiet times, these currents are usually 
limited to the daylight hours. Indeed, MAGSAT's orbit was 
located in the dawn-dusk sector in order to minimize the 
contribution of ionospheric currents to magnetic field 
observations. However, several workers have detected the effects 
on low latitude (±20°) MAGSAT data from what they believe to be 
quiet-time ionospheric currents (e.g. M«eda £JL Ai*, 1982). The 
effects have been largely on the D component in the dusk sector 
and the source has been inferred to be a meridion&l current 
system restricted to a latitude range of ±20°, and contributing a 
maximum D field variation of ±20 ranoteslas. 

The question we address here is whether appreciable 
ionospheric currc' ts develop during magnetically disturbed 
periods, i.e. during major magnetic storms. The answer is 
essential to whether one can readily use satellite data to study 
electromagnetic induction in the solid earth. The satellite 
trajectory is well above any ionospheric current system.. Any 
attempt to downward continue satellite field data through the 
ionosphere will be biased by the effects of the current sheet 
source. If the current system is known and constrained, one m<:.y 
be able to compensate for this bias. However, if the current 
system is unknown, then satellite induction studies may not be 
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feasible because of the presence of this unknown bias. 

The presence of ionospheric currents can in principle be 
detected by comparing ground-based and satelli.e measurements (a 
quantitative basis for this approach is presented in Hermance, 
1982 and Hermance^ 1983). A simple geometric model which 
illustrates our approach is presented in Figure 15. A time- 
varying sheet current source in the ionosphere induces oppositely 
directed image currents in the earth. At the surface of the 
earth a simple application of the right-hand rule indicates that 
the horizontal field from each current system (the source and the 
induced currents) will combine constructively, which results in a 
total horizontal field which is generally from I.'' to 2.0 times 
the horizontal source field alone. 

However at satellite altitudes a similar application of the 
right hand rule will show that the . *o horizontal field 
contributions combine destructively and tend to cancel. In this 
case the total horizontal field at satellite altitudes is as 
little as 5% to 20% of the horizontal source field alone. 

Hence, comparing ground and satellite observations of the 
horizontal field is a sensitive means to detect ionospheric 
current systems. Returning to our inspection on Figures 13 and 
14 we see little evidence for the presence of systematic offsets 
in either Z or H for any length of time between our satellite and 
ground-based data. In particular we feel that these data suggest 
that during the recovery phase of the storm (and perhaps during 
the main phase as well) there is little evidence for significant 
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ionospheric current flow on the disturbed Z and H components (our 
analysis of D is currently in progress). This implies that 
large-scale zonal ionospheric current systems are minimal during 
these times. Our pending analysis of the D component will allow 
us to evaluate, in a future paper, the effects of meridional 
current systems during major magnetic disturbances. At present, 
however, we feel that our analysis to-date is encouraging 
regarding the use of satellite data to study electromagnetic 
induction in the solid earth, particularly those studies which 
rely primarily on observations of the Z and H disturbed 
components . 
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STORM-TIME ELECTROMAGNETIC INDUCTION EFFECTS 
Primitive ElectrQma 9 netic Response Functions 

Data from magnetic storms during the MAGSAT mission having the 
quality shown in Figures 13 and 14 are limited. Although a 
number of magnetically disturbed periods were experienced, few 
major magnetic storms occurred. Hence at this point we feel that 
a refined analysis of the data to determine an electromagnetic 
response function is not warranted. On the other hand# an 
analysis directed at extracting some of the first order features 
of the effects of induction may be quite appropriate. To do this 
we calculate what we term to be a "primitive" response function. 
Essentially this is a gross# wide-band induction parameter 
determined from the recovery phase of the magnetic storm# where 
short term temporal fluctuations are minimal and where the source 
field can most likely be represented by a first order surface 
harmonic. 

To assess global induction effects we assume that# to the 
first-order# the earth is spherically symmetric and is 
illuminated by the magnetic field of an equivalent ring current 
system. Therefore# the source field# is assumed to be 
uniform and polarized perpendicular to the ecliptic plane over 
the earth's diameter. We adopt the usual spherical coordinate 
conventions such that in a geocentric coordinate system is the 
geomagnetic colatitude and r is the radial distance. The effects 
of asymmetries in the ring-current and gradients in the source 
field are assumed to have second order effects on the results of 



our discussion. 
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It is well-known that under these assumptions (e.g. Chapman 
and Bartels^ 1940) the magnetic field relations can be derived | 

j 

from the gradient of a scalar potential (B « -VU ) , where 1 

i 

i 

U(r,e) = (r^/r) ^ (r/r^) } cos 0 , j 

1 

» 

and are the so-called internal and external potential field | 

coeffici respectively; r^ is the radius of the earth# r is j 

the radial position of the observer and 6 is the colatitude. 

For the purpose of our present discussion we will approximate 
a region of high conductivity in the earth's interior as a super- 
conductor (o = ®) having a radius a. The magnetic field 
component normal to the surface of our super-conductor, 

B^ e “3U/3r, is zero at r ■= a. 

This condition leads to 

h " ' 

o 

where = -b^. Differentiating the above expression for the 
scalar potential with respect to r and 9 , respectively, we obtain 
the conventional geomagnetic vector components 
Z = “ [1 - ^ ' 

H = [1 + 0.5(-) ^]B Sin e . 

I r o 

These expressions can be written in terms of latitude invariant 
amplitude coefficients, analogous to those used in our previous 
discussion, above, such that 
Z = Z cos 6 , 

H = H sin 6 
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The Z and H terms rjrrespond to the amplitude coefficients# 
introduced earlier. Upon taking the ratio of Z to H we obtain 
the electromagnetic response parameter 

Z/H = - 1^- 

[1 + 0.5(a/r)-^] 

which a number of workers have employed in ground-based magnetic 
studies (e.g* Banks# 1969). 

Returning now to inspecting Figures 13 and 14# we see that the 
Z and H components both show the well-known exponential decrease 
in amplitude over a number of days during the recovery phase of 
the storm. In addition there are some shorter term fluctuations 
during which Z and H track quite well# with little phase shift. 
The close visual correlation between the Z and H time series 
prompts us to plot Z and H as a correlogram# or scatter diagram. 

A linear trend for the entire storm was subtracted from, each data 
set# using a linear least squares fit to the first 12 hours and 
the last 12 hours of the satellite data shown in Figures 13 and 
14. A correlogram of the detrended simultaneous Z(t^) and H(t^), 
sampled at various times# t^, for the recovery phase of the storm 
is shown in Figure 16. A linear relation between Z and H is 
evident# which is to be expected. The slope of the best-fitting 
straight line through these points# <Z/H>, is precisely the 
primitive electromagnetic response parameter which we will employ 
in our subsequent analysis. We also note little systematic 
difference between dawn and dusk passes. 
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These data are further segregated in Figure 17, where dawn and 
dusk passes are subdivided into passes primarily over the ocean 
(i.e. less than 10% of its path in a particular sector between 
±45^ geomagnetic latitude was over a continental land mass, see 
Figure 18) and all other passes. Again at this level of our 
analysis, there appears to be little systematic difference 
between our data subsets. 

h Simple Statistical Appioach 

By performing a mean least squares linear regression of Z(t^) 
on H(tj), where t^ represents the set of observation times, v;e 
can determine the ratio <2/H> and estimate certain statistics, 
such as the correlation coefficient. The parameter <Z/H>, our 
primitive response function, will be used to infer some of the 
gross effects from electromagnetic induction in the earth. 

The assumptions that our source field can be represented by a 
first order surface harmonic and that aliasing effects due to 
short term fluctuations between satellite passes are negligible, 
seem to be most reasonable during the later portion of the 
recovery phase of the magnetic storm we are using. Hence we will 
begin our regression at the last point of the storm and 
systematically increase the number of samples by running backward 
in time, all the while evaluating the statistical parameters 
associated with our estimate of <Z/H>. 

The late phase of the storm in which we compare samples 
reversed in time suggests a distinct correlation of two magnetic 
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channels H(tj^) and Z(tj^) which we denote# respectively# by the 
usual convention and y^. We assume that data points in this 
phase are correlated with coefficient # and have a bivariate 
normal distribution whose probability density function# for all 
physically possible storms# is given by: 


P(x,y) 


exp{- [ (x/a„) (2pxy/o..a„) + (y/a„) ^ ]/ [2 (1 -pM 1 } , 

2iTa a /I^ ^ y y 

X y 


where o and a are the standard deviations of x and y 

* y 

observations# the mean assumed to be 0. The correlation estimate 


t=t 
( I 
t=t 


n 


Vt> 


p = 


Zx^Zy2 


is a statistic whose Fisher 


2 -transform 


z = in 1 (i±2i . 

2 (1-0) (1+p) 

is a statistic with a standard normal distribution of mean = 

In and variance = 1. 

2 i— P 

In time# the statistic z tends to increase to a maximum as 
more points are included in averaging the various estimates 
toward the beginning of the storm. Qualitatively# this implies 
that the significance of the correlation increases up to this 
point. That is# the probability that the observed points could 
be produced by a weakly correlated or uncorrelated signal pair is 
decreasing to this point. 
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Howeverr as data from earlier phases of the storm are included 
in the averagesr a decrease in correlation occurs^ leading us to 
conclude that the decline in the 2 -statistic may be implying the 
modification of our basic model - perhaps the configuration of 
the source field is more complicated and our simple model is no 
longer applicable. Therefore we could reject any points after 
this local maximum, even though the correlation begins to 
eventually rise again as different phenomena become dominant. 

Given the hypothesis that p ~ p observed at the cutoff, we may 
generate confidence intervals for p which will in turn give a 
(crude) error bar for the slope of the correlation, by 
considering the values of the slope corresponding to cutoff 
points within the confidence interval for p . 

The results of this simple analysis are presented in Table 2 
for various sets of satellite passes. Our primitive induction 
parameter, <Z/H>, is the slope of the best-fitting line through 
the points shown in the scatter diagrams of Figures 16 and 17. 

The data analysed spans the period from 6 hours after the main 
phase maximum throughout the recovery phase of the storm. For 
comparison in Table 2 we indicate the actual values of the 
magnetic field amplitudes at this time (approximately 6 hours 
after the main phase maximum) , the ratio of the mean amplitudes, 
Z/H, the primitive induction parameter <Z/H> with the estimated 
range of uncertainty, and the respective correlation 


coef f icients. 
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Ui£ Import anc£ mi source field Effects opposed to Lateral 
Conductivity Variations 


In studying Table 2 , one obtains the impression that the 
difference between dawn and dusk estimates of <Z/H> (.37 ± .02 
versus .42 ± .02) seem to be significantly larger than the 
variation between passes over open ocean and passes over other 
regions. At this point therefore we are led to conclude that 
source field effects, e.g. the well-known dawn-dusk asymmetry in 
magnetic field amplitudes and perhaps other source field 
configurations, is a major contributor to the uncertainty in 
obtaining valid estimates of <Z/H>. 

A Globally Averaged Conductivity Model 


Clearly, the uncertainty in our induction parameter does net 
warrant a refined interpretation, but it may be appropriate to 
make some general observations. For example, what does a value 
of approximately 0.40 (.37 to .42) for <Z/H> imply about the 
actual conductivity structure of the earth? 


Upon rearranging our theoretical relation between Z/H and 
(a/r) given above, we obtain 


a/r 


r l+ (Z/H) 1 

^1 - 0.5 (Z/h) ^ 


We may now substitute our estimate <Z/H> for the Z/H term in this 
expression and solve for a/r. Note that the value of a/r is 
relatively insensitive to uncertainties in <.Z/H>, because of the 
cube root. In other words, the radius of an equivalent 
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superconductor is a robust parameter to estimate from our data 
base. 

The result of such calculations are summarized in Table 3 for 
dawn and dusk passes, and for the averge of both. The 
appropriate r, of course, is the radial position of the satellite 
at a nominal altitude of 400 km, so that r = 6770 km. 

Surprisingly, in view of the grossly simplifying assumptions 
we have made, the radius of an equivalent superconductor seems to 
be well estimated at a value of 5,370 (i 120) km, corresponding 
to a depth below the earth's surface of 990 (± 120) km. 
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A SEARCH FOR OCEAN EDGE-EFFECTS ON SATELLITE INDUCTION STUDIES 
Zll£ Qcgan Coast Effect 

One of the more profound examples of lateral heterogeneities 
is the difference in conductivity between oceans and continents. 
This sharp boundary leads to intense concentrations of electric 
currents, induced by transient magnetic field variations. The 
resulting secondary magnetic fields give rise to what is 
classically known as the "coast" or "ocean" effect by ground- 
based electromagnetic workers. We've chosen the southern shore 
of Australia as a region in which to observe this phenomenon and 
the storm-time enhanced westward directed ring current as our 
external inducing source. Thus we have a combination of a low 
order inducing field, a region where unhindered east-west global 
current circulation is possible, and a large scale perturbing 
structure of relatively simple geometry contained wholly within 
this region of el ectromagnetically induc.d current flow. It is 
an ideal setting in which to observe the coast effect. 

Figure 19 is an example of MAGSAT total field magnitude data 
from two contiguous passes over central Australia (passes are 
considered contiguous if they are within two degrees longitude of 
one another). In the bottom portion of the figure, one of the 
satellite tracks is plotted against a backdrop of a rotated 
Mercator map projection of the world. It serves convenient 
reference . 


An order 14 standard reference field has been subtracted from 
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the data. In addition, we've removed spikes which are apparently 
due to satellite orientation. The data at the bottom were taken 
during a magnetic quiet period selected on the basis of Kp and 
Dst indices. On the top we have data obtained during the main 
phase of a major storm during the NAGSAT epoch, beginning on 
February 14, 1980. The highly disturbed character of the top 
trace compared to the quiet time data is quite evident. In 
particular, we notice that in mid and lower latitudes, the field 
has a smooth well-behaved low order shape that one would expect 
from a ring current source. We also note that a well-defined 
dawn-dusk asymmetry is in evidence during this period of storm 
time. 

Upon closely inspecting the quiet time pass, we can see small 
scale perturbations throughout the record They range from 5 nT 
to 20 nT in magnitude. In many cases, they are due to static 
magnetization of the lithosphere. We feel the resolution of the 
MAGSAT data is consistently of sufficient quality to show 
induction effects of these magnitudes. Hence the lithospheric 
anomalies must be removed, They are separable from among the 
perturbations by meany of. two ciriteria: 

1) They must be repeatable from pass to pass ever the same 
region. 

2) leing insensitive to external time varying fields, their 
magnitudes must be invariant between contiguous quiet and 
disturbed passes. 


We turn now to the vertical field (Z) data, the Z component 
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being particularly sensitive to the ocean edge-effect. In 
Fig. 20, we have plotted this component from five contiguous 
quiet time passes over Australia. The scale is enlarged about 
six fold beyond the scale in the previous figure. We note 
remarkable repeatability of the data in the mid latitudes, and to 
some extent even in the auroral regions. Correlation of small 
scale features of magnitudes as low as 5 nT, such as over 
Australia and over Brazil, is evident. 

In Figure 22 , we compare a quiet time pass (bottom-trace) , a 
disturbance time pass from February 16 (middle-trace), and the 
difference between the two (disturbed minus quiet; top trace) . 

The top trace, therefore, shows the residual of the storm time Z 
data after removal of lithospheric effects via subtraction of the 
quiet time pass. 

At low and mid latitudes the field is dominated by enhanced 
ring current effects, as is clear from the smooth low order trend 
of the Z component. 

We see no indication of fluctuation in the data when the 
satellite crosses the south shore of Australia that would 
characterize the presence of an ocean edge-effect. This pass vas 
one of the moct disturbed during the entire MAGSAT mission. 

It is well established, of course, that an ocean effect does 
indeed exist near the coast of Australia. What our analysis 
presents therefore is a constraint on the scale size of this 
phenomenon. That is to spy, induced currents flowing in the 


"23- 


ocean adjacent to Australia have a characteristic spatial 
wavelength that is small enough so that the effect is attenuated 
below the noise level at satellite altitudes. We may obtain an 
estimate for the maximum wavelength of the ocean effect at the 
surface in the following way. 

First we note that the coast effect as observed at ground 
level along the south coast of Australia results in an 
enhancement of the Z component adjacent to the coast of at least 
a factor of 2 relative to observations 300 km or more in] and 
(Parkinson and Jones, 1979) . As we noted above, such an 
anomalous perturbation near the coastline in Australia is not 
apparent in the MAGSAT data in Figure 21. In fact, if a 
perturbation were present, it must have an amplitude less than 
10-20% of "normal" Z. We are seeing, therefore, a geometrical 
attenuation of the coast effect at MAGSAT altitudes by at least a 
factor of 10 to 20. If is typically the longest spati’l 
wavelength of the lateral dimensions of the coast effect, then 
one might expect an attenuation at a height, h, which behaves as 
exp{ -27Th/ig) . At a nominal satellite altitude of 400 km, this 
expression decreases by a factor of 10 to 20 for spatial 
wavelengths shorter than 1,100 km and 840 km, respectively. 
Typical scale dimensions of the c^-'st effect observed at the 
ground in land areas apparently are well below these values 
(generally less than 300 km) . Therefore, one can conclude with 
reasonable confidence, based on both the analytical expression 
indicated above, as welj. as by inspection of actual MAGSAT data 
(e.g. Figure 21), that one of the most profound electrical "edge- 
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effects" which influence ground-based magnetic induction studies 
has a minimal contribution to magnetic field components at MAGSAT 
altitudes. 
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CONCLUSIONS 

The work summarized here is best described the 
reconnaissance phase of evaluating the contribution that 
satellite observations can make to studying electromagnetic 
induction in the solid earth. Although certain elements of our 
analysis are qualitativef several points can nevertheless be 
made. First , satellite data apparently suffers far less from the 
effects of near surface lateral heterogeneities in the earth than 
do ground-based data. Second . zonal ionospheric currents during 
the recovery phase of major magnetic storms appear to be minimal, 
at least in the dawn and dusk sectors where KAGSAT was flown. 
Hence the internal contributions that satellites observe during 
these times is in fact due primarily to induction in the earth 
with little or no contribution from ionospheric currents. Third . 
the interpretation of satellite data in terms of primitive 
electromagnetic response functions suggests that more 
quantitative interpretations hold great promise for refining 
models of the electrical conductivity in the upper mantle. Even 
the grossly over-simplified interpretation presented here results 
in a surprisingly well-resolved radius for an equivalent super- 
conductor representing the conductivity region of the earth's 
interior (5,380 + 120 km). This is indeed encouraging for 
follow-on work in this area. What is required now is a refined 
spectral analysis of a number of different transient events, 
accounting for a more complicated source field than has been 
assumed hitherto. Having resolved the frequency dependence of 
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the electromagnetic response parameter r one will be able to 
investigate the global conductivity structure o£ the earth's 
interior using more realistic models than the one used here. 
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Table 1 . Calculated offsets between ground-based 
and satellite data. 


Residual = Smoothed satellite data - ground-based data 


Component 

Residual 



Dawn Sector 

Dusk Sector 

H 

-24.3(±b,0} * 

-19.8(110.6) * 

Z 

24.1i-'9.6)* 

16.9(110.4) * 


* ±1 standard de ;iation 


Tabic 2. Estimates of the primitive induction parameter 'Z/H> 
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Table 3 . Range of simple conductivity models 
for the earth based on primitive 
response function from satellite data. 



Response 

Radius 

of Equivalent 


Parameter, <Z/H> 

Super 

-Conductor 

Dawn 

-.37 

(±.02) 

5,480 ( 

±75) km 

Dusk 

-.42 

(±.02) 

5,300 ( 

±80) km 

"Average" 

-.40 

(±.03) 

5,370 ( 

±120) km 


-27- 


ORIGINAL PAGE IS 
OF POOR QUALITY 


REFERENCES 

BankSr R.J.» Geomagnetic variations and the electrical 

conductivity of the upper mantle^ Geophvs. J. Eoy. Asti on « 

Sqc. . 11, 457-487, 1969. 

Chapman, S. and J. Bartels, Geomagnetism * Clarendon Press, 

Oxford, 1940. 

Didwall, E.M., The electrical conductivity of the Earth's upper 
mantle as estimated from satellite measured magnetic field 
variations, Ph.D. thesis, Johns Hopkins University, Baltimore, 
MD, 1981. 

Didwall, E.M. The electrical conductivity of the upper mantle as 
estimated from satellite magnetic field data with implications 
for temperature. Submitted to J. Geophys. £££., 1983. 

Hermance, J.F., Electromagnetic induction and source field 
effects in MAGSAT data (Abstract) , Programme aM Ab.S.tXACL £ , 
Fourth lAGA Scientific Assembly, edited by N. Fukushima, p. 
538, lUGG Publ. Office, Paris, France, 1981a. 

Hermance, J.F., Are there induction effects in Dst corrections to 
MAGSAT data? (Abstract), £QS, £2, P- 269, 1981b. 

Hermance, J.F., Finite source fields coupled to lateral 

conductivity heterogeneities; Effects on magnetotelluric and 
magnetic gradiometric deep sounding experiments (Abstract) , 
Programme and Abstracts * Fourth lAGA Sciei.tific Assembly r 
edited by N. Fukushima, p. 566, lUGG Publ. Office, Paris, 
France, 1981c.. 

Hermance, J.F., Model simulations of possible electromagnetic 

induction effects at MAGSAT altitudes, Geophys . Res . LeUt*, £, 
373-376, 1982. 

Hermance, J. F. and M. Rossen, Regional and global 

electromagnetic induction effects in MAGSAT satellite data 
(Abstract) , Sixth Workshop on Electromagnetic Induction in the 
Earth and Moon, Working Group I-S, lAGA, University of 
Victoria, British Columbia, Canada, 1982. 

Hermance, J.F., Electromagnetic induction by finite wave-number 
source fields in 2-D lateral heterogeneities; The transverse 
electric mode. Submitted to Geophys . J.R. astr . 1983. 

Langel , R.A., Internal and external storm-time magnetic fields 
from spacecraft data, I AG A Program Abstracts loi. XYI Genei.al 
Assembly . Genoble, France, 1975. 

Langel, R.A. and R.H. Estes, Large-scale, near-earth, magnetic 
fields from external sources and the corresponding induced 
internal field, NASA Technical Memorandum 85012, Goddard Space 



-28- 


Flight Centeif Aprils 1983. 

Langelf R.A.f R.H. EsteSf G.D. Mead# E.B. Fabianor and E.R. 

Lancaster r Initial geomagnetic field model from MAGSAT vector 
data, Geophys . Res . Lett . » 2, 793-796, 1980. 

Langel, R., J. Berbert, T. Jennings, and R. Horner, MAG5AT 

processing ; A report for investigators . Technical Memorandum 
82610, 178 pp., NASA, 1981. 

Langel, R.A., The magnetic earth as seen from MAGSAT, initial 
results, Geophys . £££. Lett . . £, 239-242, 1982. 

Langel, R.A. , G. Ousley, J. Berbert, J. Murphy and M. Settle, The 
MAGSAT mission, Geophys . Res . Lett . , 2, 243-245, 1982. 

Langel, R.A., J.D. Phillips, and R.J. Horner, Initial scalar 
magnetic anomaly map from MAGSAT, Geophys . Res . Lett . , 2, 
269-272, 1982. 

Langel, R.A., C.C. Schnetzler, J.D. Phillips, and R.J. Horner, 
Initial vector magnetic anomaly map from MAGSAT, Geophys . £££. 
Lett , a, 273-298, 1982. 

Maeda, H,, T. lyemori, T. Araki and T. Kamei, New evidence of a 
meridional current system in the equatorial ionosphere, 
Geophys . Res . Lett . , £, 337-340, 1982. 

Parkinson, W.D. and F.W. Jones, The geomagnetic coast effect. 

Rev . Geophys . And Space £hys.» 12f 1999-2015, 1979. 


-29- 


ORIGINAL PAGE IS 
OF POOR QUALITY 


FIGURE CAPTIONS 

Figure 1. Locations of standard magnetic observatories used in 
this study. 

Figure 2. Ground-based vector component magnetic data from the 
Tucson standard observatory. Hourly averages were low pass 
(Hann) filtered such that full amplitudes are present at a 
12 hour periods the amplitudes are reduced to one-half at a 
6 hour period and to zero at a 3 hour period. The arrow 
labeled "NDNT" denotes the local midnight values which were 
used to establish the zero baseline. 

Figure 3. Reduced global ground-based data for the H-component 
during the dusk, midnight, and dawn sector. The time 
extends from 00:00 hrs. Feb 12 to 00:00 hrs. Feb 23, 1980. 

Figure 4. Reduced ground-based data for the Z and H components 
during the dawn sector, 00:00 hrs. Feb 12 to 00:00 hrs Feb 
23, 1980. 

Figure 5. Geographic location of satellite passes used in this 
study. Geomagnetic latitudes of ±45^ are shown for 
comparison. Data gaps are shown as interruptions in each 
satellite track. 

Figure 6. H-field amplitude coefficients for satellite data 

(solid pts.) compared with ground-based data (open circles). 
Note the similarity in form, but the presence of a small, 
persistent offset. 

Figure 7. Z-field amplitude coefficients for satellite data 

(solid pts.) compared with ground-based data (open circles). 

Figure 8. Histograms of offset residuals during the dusk sector 
for the H component satellite data minus the ground data. 
Top: two days at the beginning; Middle: two and a half 

days at the end; Bottom: the total of both beginning and 
end. 

Figure 9. Histograms of offset residuals during the dawn sector 
for the H component. See caption for Figure 8. 

Figure 10. Histograms of offset residuals during the dusk sector 
for the Z component. See caption for Figure 8. 

Figure 11. Histograms of offset residuals during the dawn sector 
for the Z component. See caption for Figure 8. 

Figure 12. Z-component amplitude coefficients (top) and H- 

component amplitude coefficients (bottom) for the satellite 
data during the dawn secroi with constant offset removed. 

Figure 13. Comparison of fully reduced satellite data (solid 
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points) with ground-based data (open circles) for the dawn 
sector. 

Figure 14. Comparison of fully reduced satellite data (solid 
points) with ground-based data (open circles) for the dusk 
sector. 

Figure 15. The geometry of external ionospheric source currents 
and internal induced currents. 

Figure 16. Correlogram of simultaneous Z and H amplitudes for 
all satellite passes during the recovery phase of the storm 
from approximately 6 hours after the main phase maximum to 
00:00 hr. Feb 22 , 1980. Results for the dusk sector are 
shown as solid points, those for the dawn sector are open 
circles. 

Figure 17. Segregation of Z and H amplitudes into correlograms 
for passes over open ocean (OCN) and all other pases (ALL) . 
These are shown for the dusk and dawn sectors, to see if a 
source effect is evident. 

Figure 18. Those satellite tracks classified as ocean passes 
between geomagnetic latitudes of ±45^. 

Figure 19. MAGSAT total field magnitude data from two contiguous 
passes (within 2^ of each other) for the satellite track 
shown. 

Figure 20. MAGSAT vertical field (Z) data from five contiguous 
quiet time passes. Scale is enlarged about 6 fold relative 
to Figure 19. 

Figure 21. Comparison of MAGSAT vertical field (Z) data for a 
quiet time pass (bottom-trace), a disturbance time pass 
(middle-trace) and the difference between the two (disturbed 
minus quiet; top-trace). 
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ABSTRACT 

Interest in electromagnetic induction in a laterally 
heterogeneous earth by a finite source field spans a variety 
of geophysical disciplines. Most efforts reported to date are 
somewhat restricted in that certain scaling relations cannot 
be violated, or else the modelling results become invalid. 

For example certain approaches are appropriate only when the 
dimensions of the heterogeneity are small compared to the scale 
size of the source field. In addition, interest has recently 
developed in using satellites to study induction in the earth. 
The particular case when the satellite is above the source 
(such as for ionospheric current systems) is one which poses 
special problems to several existing approaches. 

This paper describes a more general approach to solving 
these problems and broadens the range of scaling relations 
over w..xch such models can be applied. The specific geometry 
considered is an equivalent 2-D current source directed parallel 
to the strike of a 2-D anomalous structure within the earth 
having a characteristic dimension of L; this geometry rep- 
resents the transverse electric (TE, or E-parallel) induction 
mode. 

A finite source field, which is located at some position 
within the modelling region, is simulated by discontinuities 
in the magnetic (or electric) field appropriate to the current 
flowing in the source. However, the fields reflected from the 
earth are continuous through the source region. Such an 
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approach not only has certain mathematical advantages, but 
also allows us to evaluate the electromagnetic fields above 
current sources in the ionosphere such as might be observed 
by satellites. 

The case when the source field itself has a finite dimen- 
sion ir.poses additional considerations for the problem of 
designing adequate 2-D models. The mesh spacing must be 
sufficiently dense to allow quadratic forms employed by the 
numerical algorithm to adequately represent trigonometric 
functions in the horizontal direction and exponentially damped 
functions in the vertical dimension. Moreover, in some cases, 

close to electrical discontinuities the induced field may behave 

2 

more like 1/r or 1/r . This often necessitates excessively 
close nodal spacing near electrical discontinuities. Whereas 
for the case of uniform source fields incident on a 2-D earth, 
one often increases the nodal spacing to very large values as 
one approaches the edges of the model, this is generally in- 
appropriate for finite source fields. 

Our results suggest that while fields may be perturbed 
in the close vicinity of lateral heterogeneities, at even 
modest offset distances from the contact the fields are more 
influenced by what is beneath the observer than by what is at 
some lateral distance away. This implies that in many cases, 
simple 1-D plane-layered models may be appropriately used in 
interpreting field data in a "piece-wise" sense, even when 
the total structure is more complicated laterally. 
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The fact that the fields induced by finite sources tend 
to decouple fairly rapidly with distance to either side of a 
vertical discontinuity prompted us to consider the effects that 
this phenomenon has on magnetic gradiometric experiments. Our 
results are encouraging in that, at least in some cases, regional 
gradients of the horizontal field (even where the sites span a 
major discontinuity) can be combined with local values of the 
vertical field, to yield reliable estimates of electromagnetic 
response functions. 

At satellite altitudes, the difference between the results 
for the uniform source field and those for a finite source 
field in the ionosphere are profound. The total horizontal 
field amplitude produced by the uniform source is approximately 
unity, suggesting that for remote sources in the magnetosphere 
the horizontal field is insensitive to the effects of lateral 
variations in the earth's conductivity. On the other hand, 
the total horizontal field amplitude produced by the finite source 
field in the ionosphere is reduced by more than an order of mag- 
nitude. This is due to the fact that the satellite is above the 
source so that the induced and the source horizontal field com- 
ponents tend to cancel and the resultant H-field is quite small. 

In contrast, the finite source produces an anomalous vertical 
field at satellite altitudes which is of the same order of mag- 
nitude as that produced by the uniform source. 
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INTRODUCTION 

Motivation 

This paper considers electromagnetic induction in a 
laterally heterogeneous earth by a source field having finite 
dimensions. Interest in this problem spans a variety of geophys- 
ical disciplines which range from exploration studies, using 
artificial source fields, to global induction studies, using 
natural sources in the ionosphere and magnetosphere. 

While a number of investigations have been reported in the 
literature which deal with various aspects of this problem, most 
efforts are somewhat restricted in that certain scaling relations 
cannot be violated, or else the modelling results become in- 
valid. For example certain approaches are appropriate only 
when the dimensions of the heterogeneity are small compared to 
the scale size of the source field. 

In addition, interest has recently developed in using satel- 
lites to study the coupling of finite source fields to the con- 
ducting earth (e.g., Hermance, 1982). The particular case when 
the satellite is above the source (such as for ionospheric current 
systems) is one which poses special problems to several existing 
approaches and has not been studied up to the present time ex- 
cept for simple, plane-layered earth models. 

The following discussion describes our progress in developing 
a more general approach to solving these problems. We not only 
want to broaden the range of scaling relations over which our 
model can be applied, but we are especially interested in ad- 
dressing those problems related to the effects of induction as 


- 2 - 


seen at satellite altitudes (Hernance, 1982). The geometry 
considered in this paper is illustrated in Figure 1, which 
shows an equivalent 2-D current source directed along the 
y axis and flowing in the ionosphere at some height above the 
earth's surface. The source current varies in time according 
to exp(ifa)t) and flows parallel to the strike of a 2>D anomalous 
structure within the earth having a characteristic dimension 
of L. In the parlance of theoretical electromagnetism, this 
geometry represents the transverse electric (TE, or E-parallel) 
induction mode. 

The procedure we follow involves establishing analytical 
relations between field quantities at the boundaries of our 
modelling region and then solving for the fields within our 
modelling region using a numerical algorithm. In principle 

this is a standard approach to such problems. 

Previous Work 

In the past, some workers have preferred an algorithm 
whereby the uppermost boundary of the modelling region (shown 
as the dashed line in Figure 2) is actual... below the source, 
and the value of the field is constrained to represent the sum 
of that of the incident source from above and that of the in- 
duced field from below. To account for the induced field one 
needs to introduce an a priori value for the average reflectance 
of the earth. This in general is most appropriate when the 
characteristic spatial wavelength of the source (X_) as well as 
the height of the top-most boundary is much greater than the 
characteristic dimension (L) of the anomalous region. 
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Other workers have preferred boundary conditions such 
that in the vicinity of the source, one **sees'* the effects of 
the source field alone (Figure 3), while at "infinity”, the 
electromagnetic fields go to 0. These conditions are appro- 
priate when Xg is small compared to the characteristic dimension 
of the model - in some sense, this is the opposite extreme of 
the case illustrated in Figure 2. 

Swift (1971) used a transmission surface analog of a 
finite difference approach to simulate induction by a line 
current in an imbedded heterogeneity. He imposed a plane-' 
wave characteristic impedance boundary condition of 

* (iup/o)^^^ 

at the top, sides and bottom of the model. Although this ap- 
proximation is correct for some cases, it is not correct in 

general. The impedance (defined as -E /H ) of a homogeneous 

7 ^ 

earth to a line current source flowing in the y-direction is 
generally not the same as the impedance of a homogeneous earth 
to a plane wave source (see for example. Wait, 1953). In the 
case where the impedance of the earth to either source is ap- 
proximately the same (i.e. when 6<< Xg, 6 being the depth 
of penetration) , then these boundary conditions are appro- 
priate when the distance to the boundaries is large when com- 
pared with the characteristic dimension of the source field as 
well as when compared with the characteristic dimension of the 
lateral heterogeneity (or the depth of penetration). 
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Stoyer and Greenfield (1976) also used plane-wave char- 
acteristic impedance boundary conditions at the edges of their 
model. They found that their results within the interior of the 
model, close to the source, appeared to be reasonably reliable 
when compared with analytical (but numerically integrated) 
Solutions. However when the edge of the model was approached, 
the quality of their results degraded appreciably. 

Coggon (1971; p. 140) felt that when simulating induction 
by a line source parallel to a two-dimensional heterogeneity 
(E-polarization), the appropriate condition to apply at ex- 
terior boundaries is that the secondary electric field is 
zero. This, of course, also implicitly assumes that either 
the scale size of the source is small or that the dimensions 
of the lateral heterogeneity are small. 

Hibbs and Jones (1973) adapted the numerical algorithm 
of Jones and Pascoe (1971) and Pascoe and Jones (1972) to 
treat the case of induction by finite source fields in a two- 
dimensional current system flowing parallel to the strike of 
conductive heterogeneities (E-polarization). An implicit re- 
striction on their model, although one they did not emphasize, 
was that the plane-layered structure must be the same at either 
vertical edge of their exterior boundaries. Moreover, in the 
example they discussed, they assumed uniform reflectivities at 
an altitude of only 50 km which implies the effects of lateral 
heterogeneities are minimal at this height. This assumption 
is appropriate if the scale size (L) of the lateral heterogeneity 
is small compared to the height at which this condition is 
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applied. In the case they considered, the lateral heterogeneity 
must have a dimension L << 50 km. 

Quon e^ a^, (1979) have considered the effects of what are 
probably extreme cases of dipole and line current sources at 
heights of only 110 km. Although they do not discuss their 
numerical algorithm in detail, they describe an example of 
induction by finite sources in an earth which consists of 
two quarter-spaces (of resistivity 100 and 10 ohm-m, respec- 
tively) separated by a vertical contact. Each sub-region 
extends to infinity in its respective quarter- space. 

Weaver and Thompson (1972) have shown for a restricted 
case the importance of properly selecting the dimensions of 
a numerical model and specifying the boundary conditions at 
its edges. They point out that improperly posed boundary 
constraints often perturb numerical results over significant 
regions of the model's interior. 
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Present Model 

The case to be considered here is somewhat more generalized 
than much of the previous work. As illustrated in Figure 4, 
we want to simulate a finice source field which is located 
at some position within the modelling region. The presence 
of the source is simulated by discontinuities in the magnetic 
(or electric) field appropriate to the current flowing in 
the source, but the fields reflected from the earth are con- 
tinuous through the source region. Such an approach has certain 
mathematical advantages, but also allows us to evaluate the 
electromagnetic fields above current sources in the ionosphere 
such as might be observed by satellites. Alternative formula- 
tions, such as the class of models represented in Figure 2, are 
sometimes not compatible with this application. 

At the top-most and bottom-most boundaries, characteristic 
impedance conditions are applied (accounting, of course, for 
the finite wave-number of the source). In addition, the earth 
is assumed to be plane layered at either of the side boundaries 
(but the vertical structure can be different at either side). 

In the interior of the model, the theory allows for quite 
arbitrary 2-D lateral heterogeneities. Clearly, however, it 
is required by the physics of the problem that all boundaries 
be far from the localized effects of fields associated with 
electrical discontinuities within the model's interior. In 
the following sections we consider in detail the implementation 
of these concepts, first in an analytical sense and then in 
terms of numerical operators. 
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BOUNDARY CONDITIONS AT THE 
LEFT AND RIGHT EDGES OF THE MESH 

Plane-Layered Earth 

For the class of problems we wish to consider in the pre- 
.sent analysis, we follow the conventional assumptions employed 
for two-dimensional models using unif or n- source fields (Madden 
and Swift, 1969; Wright, 1970; Jones and Price, 1970): the 

ear h is plane-layered at the extreme left and right vertical 
edges of the mesh (i.e. the outer-most column of nodes); and 
that these edges are sufficiently far from scattering centers 
that one may neglect the effects of lateral heterogeneities 
on the fields at the edges. We wish to assume, however, that in 
general the horizontal layers along the left edge may be dif- 
ferent from those along the right edges. This assumption allows 
us to apply Dirichlet boundary conditions at the edges approximate 
for a finite source above a plane layered earth, a problem which 
is solvable analytically (Price, 1950; Wait, 1953; Hermance and 
Peltier, 1970; Peltier and Hermance, 1971, Hermance, 1978). 

We now describe the application of these boundary conditions 
to our finite difference simulation. 
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A Surface Impedance Recursion Relation 

Expressions for determining the dovmward- looking impedance 
at an arbitrary surface above a n- layered earth have been 
derived by many workers (e.g. Wait, 1953) and the problem need 
not be discussed in detail here. However for the sake of con- 
sistency in notation etc., we briefly outline the development. 

To derive the impedance recursion relation, we begin by 
considering the orthogonal, electric and magnetic field com- 
ponents, and Hy, and assume E^ and E^ are 0. For self- 
consistency, the electric and magnetic fields must be related 
through Maxwell’s equation, for example 

VxE = -3B/at , (1) 

or 

(VxE)y = -iupHy , C2) 

where 

(VxE)y « . (3) 

Therefore 

Hy « -(l/iwp)3E3^/3z . (4) 

It', in the jth medium, the electric field follows a relation of 
the form 

Exj • + Bje‘^j^)cos6x , (5) 

1/2 

where 0 = k(x-x ) and y. * (k^ + iwpa- - w^ye) ' , then 

A U J J 

Hyj * -(Yj/iuy) [A^e^j • B^e ]cos6x . (6) 
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At the jth interface, we define the surface impedance, 

. (7) 

which of course is not to be confused with the generalized 
charracteristic impedance, * (ioiv/Vj)* of a strictly homo- 
geneous medium. 

'dative to a coordinate system located at the jth inter- 
face, the ration f^ormed from (5) and (6), and 

evaluated at z « 0 is given by 


f^xj-l^yj-1^ 


z®0 




( 8 ) 


which we term our generalized surface impedance at the top of 
the jth medium; in other words, denoting this parameter with a 
subscript j , we obtain 


which can be rearranged to the form 




(9) 


( 10 ) 


Expressions such as (5) and (6) can be used to express the 
generalized surface impedance at the top of the j-1 layer, 
relative to a coordinate system at the jth interface: 


j-1 




I-Adj.l 


-Z. ,{ 





i-l' i 

3 Aj.je 3-1-B^.i 


■j-1 


( 11 ) 
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where the argument of the exponent is given by 
However a relation between Aj.j and Sj.j is expressed in (10), 
so, for example, upon solving (10) for substituting the 

resulting expression into (11), rearranging terms and cancelling 
the B- . terms which appear in the numerator and denominator, 

J ‘ * 

we obtain 


(zj - Z^.;^)e‘®j-1 - (zj + 


( 12 ) 


Combining like terms in T? and 2j_^, respectively, we obtain 


. , Z^cosho. 1 + Z. .sinha- , 

zi-i = ^ ^ 

■' Z^ sinhoj _ j + Zj _ jCosha^ _ j 


(13) 


or upon dividing through by cosh a 


3 


- 1 * 



Z^ + Z. -.tanh a- , 
Z._, ^ ^ 

^ Z^tanh 


we arrive at the form 

(14) 


for the surface impedance at the top of the j-1 layer, where 
Z^ is the surface impedance at the top of the jth layer (or the 
bottom of the j-1 layer). 

Clearly we may use expression (14) as a recursion relation. 
Assuming that a homogeneous medium is present beneath the bottom- 
most interface (j = n) , we have Z^=Z^^. Using (14) we can prop- 
agate this relation to the next higher interface (j « n-1) to 
determine z'^’^. This recursive operation can be continued to 
progressively shallower depths until we reach the node directly 
beneath our source current. 
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ELECTROMAGNETIC FIELDS OF THE SOURCE 
Source Current Density 

It is well-known (e.g. Hermance, 1978) that an arbitrary 
two-dimensional source flowing in the y-direction, having a 
current density JyCx»z»t), can be expressed as a superposition 
of elementary sources having a spatial wavelength X , or wave- 
number k« Zv/X^. For the case of a thin horizontal current 
sheet flowing at a fixed height (2*2^) above the earth’s sur- 
face (Figure 5) we can expiess the current density as 

Jy(k) « I^^(k)6(2-z^)cos[k(x-XQ)l (15) 

where 6(2-2^) is the Dirac delta function and I^(k) is the 
magnitude of the coefficient associated with the kth spatial 
wave number (which we assume has the dimensions of a sheet current 
density, i.e. amperes/m). For the present discussion we will 
only consider sources of the above form, recognizing that if 
needbe, more complicated sources can be easily synthesized from 
the resulting expressions (Hermance, 1978). In the present 
context we emphasize that sheet currents such as (15), which 
may be reconstructed from ground-based magnetic field measure- 
ments, may not represent actual current systems in nature. In 
some cases they should be looked upon as only hypothetical 
current systems which produce magnetic fields at the earth's 
surface which are equivalent to the magnetic fields produced 
by actual ionospheric or magnetospheric current systems - the 
actual current systems may have quite a different geometry 


in three dimensions. As an example, even certain classes of 
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electro jets, which may in face be terminated at either end by 
field-aligned currents, generate electromagnetic field varia- 
tions over restricted regions at the earth's surface which can 
be represented as the fields of equivalent current sheets such 
as (15), Above the source of course, or within the «^ource, 

*the expressions will be completely different. However the con- 
cept of equivalent sheet current sources is a useful artifice 
for evaluating electromagnetic perturbations at the earth's 
surface when the sources are some distance above; particularly 
when the sources become quasi two-dimensional as they seem to 
do for many classes of electromagnetic phenomena ranging from 
pulsations to magnetic substorms. In order to emphasize that 
it is the effects of the current system that we are attempting 
to simulate over a limited region of interest, rather than 
attempting to physically model the '.ources themselves, the 
sheet current represented in (15) will be termed an equivalent 
sheet current density (e.g. used by Price, 1962). 
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Electromagnetic Field Components 

Directly beneath the source represented by (IS) and im- 
mediately adjacent to it, we must have, for the source field 
alone , 

» Jy/2 , (16) 

where, for convenience in this and all following equations, we 
suppress the k index whenever possible without leading to con- 
fusion in the resulting terms. 

An appropriate form for the horizontal component of the 
magnetic field observed at a position (x,z) beneath a source 
at (Xq>Zq) in a homogeneous medium is, according to Hermance 
(1978), 


exp[-Y(z-z^)]cosk(x-XQ) , (17) 

where y* (k* ••• iuyo - For convenience we will use the 

following constants to represent the arguments of the terms in 
(17): 

' yCz-Zq) , (18a) 

®x ‘ 

For the class of problems considered here, we will assume 
B * y© ii everywhere, wh®re y^ * 4irxl0 henries/m is the 
permeability of free space. In this case 

V • H ' 0 (19) 

everywhere, and we are led to the following condition on the 
derivatives: 
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aH^/dz • -3Hj^/ax . (20) 

Upon substituting (17) for and taking the derivative 
with respect to x on the right hand side, we have 

aH^/3z « (Ij,k/2) exp(- 62 )sinejj , (21) 

which can be easily integrated to the form 

“z * -fIo^/t2y)] exp(-e^)sine^ . (22) 

The electric field has only a y-component and is related to 
the magnetic field components through Faraday's Law of Induction. 
In other words, if we assume the following form for the electric 
field of the source: 

Ey ' exp(-e^)cos6^ , (23) 

the coefficient Bf is determined from 

o 

H - -(l/iwp)7xE . (24) 

Upon taking the curl of the electric field on the right hand 
side of (24), we obtain 


quality 


= (l/iwy) 3Ey/3z 
« -(Y/i(oy)B® exp(-02)cose^ , 

= (k/iuy)B® exp(-02)sine^ . 


(253) 

(25b) 

(26) 
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By comparison of (25) with (17), or (26) with (22) it is clear 
that 

- -[iuyIj,/(2Y)l . (27) 

To summarize, we assume an equivalent sheet current having 

% 

the spatial form 

where is the fixed amplitude of the spatial harmonic at the 
wavenumber k, and 0^«k(x-x^) so that the source is centered 
at x*x^. In this case the source-field electromagnetic field 
components at a distance z - z^ in a homogeneous medium beneath 
the source are given by 



- [iuvlQ/(2Y)] exp(-e^)cose^ , 

(29a) 

«x - 

(1^/2) exp(-e^)cose^ , 

(29b) 

Hz ' 

-[I^k/(2Y)] exp(-e^)sine^ , 

(29c) 

where 6 = yCz'Z ) . 

On the other hand, in the medium 

above the 


source we have the following expressions for the source field 
components : 


- [i(DUljj/(2Y)] exp(+e^)cos6^ , 

(30a) 

-(1^/2) exp(+02)cose^ , 

(30b) 

-[lQk/(2Y)] exp(+02)3in0^ . 

(30c) 
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Hence the electric field of the source has the same spatial 
dependence as the horizontal magnetic field component, whereas 
it has the opposite sense to the direction of current flow. 

A few moments thought and the right-hand-rule will justify the 
fact that is maximum at positions of x for which and 
are minimum. 


A Downward -Propagating Recursion Relation 

We assume that in the medium between each node along the 
left and right hand vertical boundaries we have the following 
expressions for the total (source plus induced) electromagnetic 
field components: 


Eyj - [Aj exp(y^z) + B^ exp(-YjZ))cose^ , (31a) 

Hxj * (Yj/iwv)[Aj exp(Y^z)-Bj exp(-Yj z) ]cos6^ , (31b) 
H_. - (k/iu)v)[A. exp(Y^z) + B. exp(-Y^ z) ] sine . (31c) 

j j j j *■ 


The subscript j refers to field values and material constants 
appropriate for the jth medium, defined in our convention as 
the region directly below the jth node, but above the j + 1 node. 

To be definite, we assume the equivalent sheet current source 
is in free-space above, but infinitesimally close to the j “ jg 


node. We also assume that the surface impedance 


at this 


node has already been calculated, hence we have the condition that 




( 32 ) 
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For the special case when ■ 1, we will assume the source is 
immediately above, but closely adjacent to the topmost row of 
nodes. In the latter situation we will assume that it is ap- 
propriate to apply Dirichlet boundary conditions along the 
topmost row of nodes which represent the total (source plus 
induced) field components above a plane -layered earth. For 
this to be possible, any lateral heterogeneities must be well 
within the interior of our finite difference mesh, the left 
and right hand conductivity structure must be plane -layered 
and identical, and the scale-size of the heterogeneity must 
be small compared to its distance from any boundary. Although 
somewhat restrictive, these conditions still permit asymmetric 
source fields (i.e. in general cos0^^® f cos6^^®) , and the 
lateral heterogeneity need not be symmetric. 

We assume that in the infinitesimally thin region between 
the source and the interface, the electric field can be 
represented by 


Ey = -[iwyI^/(2y^j)]{exp(-YQz) - R(k)exp(YQz)}cose^ , (33) 

where » (k^-u^ye) and R(k) is the reflection coefficient. 
The corresponding horizontal magnetic field is given by 




Upon taking the ratio of (33) to (34) at z = 0 and equating the 
result to the surface impedance -Z^ U we obtain 

I 


-Z 


3 



-(ia3y/YoHll-R]/[l+R]> 


(35) 
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which can be solved for the reflection coefficient 

R(It) - [(Zo- Zj)/(Z ♦z3)]. (36) 

j J s 

where • iwy/Y^ is the generalized characteristic impedance 
of free-space (i.e. includes the spatial wavenumber, k) . 

At the jg interface we equate the horizontal electric field 
components across the jg"! and media, as well as the horizontal 

magnetic field components. Using (31a) and (31b), this leads to 


E 


y 


3s 



(A^ + Bj)cose^ , 
(l/Z^)(Aj - B^)cos6^ , 


(37) 

(38) 


the left hand sides of the equations being given by (33) and 
(34) respectively. We might note at this point the general ap- 
plicability of relations of this form to any interface, such 
that the electric and magnetic field components at the bottom 
of any layer, j-1, is related to the coefficients in the layer, 
j, below, through the forms 




(Aj + Bj)cose^ , 


(1/Z.)(A. - BJcose , 

j j j * 


(39) 

(40) 


where the expressions are evaluated at the j interface. We can 
solve for the coefficients 


Aj = (1/2) (eJ' 


. ♦ z.hI 

J 3 


M 


j l/cosG^ 


9 


(41a) 
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Bj - (1/2 )(eJ'^|. - . (41b) 

Using the foTms (31a) and (31b) we can now determine the 
electric and magnetic field components at the bottom of the 
jth layer through (31a) and (31b) such that 

4li^l * exp(a^ + B. exp(-a^)]cose^ , (42a) 

and 

hJ - (l/Z^)[Aj exp(a^) - Bj exp(-ap]cos6j^ , (42b) 
where a. * y .d. . 

We now have the architecture for developing a recursion 
relation for the downward propagation of field values at each 
node. We begin at the interface directly below the source by 
determining the surface impedance, and from (36), the reflection 
coefficient. The reflection coefficient is substituted into 
(33) and (34) to yield E and H , respectively. These expres- 

y ^ 

sions are substituted into (41a) and (41b) , along with the 
appropriate characteristic impedance 2^ , to yield the coeffi- 
cients Aj and Bj . The c'Jcif icients are next substituted into 
(42a) and (42b) multiplied by the exponential factors which are 
indicated and result in values for the electric and magnetic 
field components at the next lower interface. We have essen- 
tially continued our fields downward from the jth interface to 
the j+1 interface. 

This process can be iterated until we reach the bottom of 
the grid. If this is done numerically, then it is usually 
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ad\antageous to evaluate the magnitude of the field components, 
normalized by the values directly beneath the source, say, 
and when the magnitudes become vanishingly s.all, to set nodal 
values at greater depth equal to zero. 


Propagating the Recursion Relation Above the Source 

In the medium above the source we must consider source 
terms which decrease with altitude. For the sake of simplicity, 
we neglect the effects or reflections above the source so that 
at the jth node we have 


Ey - - [iu)plQ/(2Y^)]{exp(aj) - R(k)exp(aj)}cos6^ , (43a) 

Hy “ (I./2) {-exp(o.) + R(k)exp(a.) }cose , (43b) 

AO J J A 

and 

* -[IoJ«/(2Yo)He3cp(aj) - R(k)exp(a^))sine^ , (43c) 

2 2 1/2 

where such that Yq “ (k - u uc) » and z is measured 

relative to a coordir'’te system having its origin (z*0) at 

the source, j*j„ (z increasing positive downward). These re- 

lations are simply used to calculate , H • and H . at each 

J ^ J 

node; j, above the source, after Zj is specified. 
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BOUNDARY CONDITIONS AT THE 
TOP -MOST AND BOTTOM- MOST ROW OF NODES 

Impedance Boundary Conditions 

We assume that all sources are contained within the region 
spanned by the model and that no reflections occur from above 
the top-most nor from below the bottom-most row of nodes. 

Hence the flow of power is everywhere outward through 
the top and bottom boundaries of our model and we are able to 
impose a generalized impedance condition on the fields crossing 
these boundaries. 

At the top-most boundary we impose the condition that 
E /H J - Zt • 

' 2*Z^ 

= iuv/Y^ » (44b) 

where = (k + iwyo^ - w pc) ; the subscript, t, denoting the 
material parameters associated with the homogeneous region at 
and beyond the top boundary. The positive impedance term in 
(44a) is compatible with the upward flux of power. 

We have of course, from Faraday's Law of Induction (1), 

that 


H - (l/iuy)3E /3z , (45) 

A / 

so that the co..dition (44b) can be posed as a condition on the 
electric field Ey and its vertical derivative such that, at 


Ey - (l/Y^)?iEy/3z « 0 . 


(46) 
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At the bottom-most boundary we have downward flowing power 
and impose the condition that 


Ey/»x 


- - 2 , 


Z=2, 


-iuy/Y^ 


(47a) 

(47b) 


which, upon invoking (45), can be expressed at 


(48) 





IMrLEMLNTATlCN OF THE THEORY 
USING FINITE UIFFERENCE SIMULATION 

A Genera 1 FJ n itc P ifference Operator 

The numerical solution of problems concerned with electro 
magnetic induction in heterogeneous structures can be accom- 
plished using a number of techniques (see, for example, the 
review by Ward M al_- IS'^4'). Of these techniques, we have 
chosen to employ the finite difference method. In brief, the 
present formulation has the following four features. 

First . complex field quantities are calculated directly. 
This facilitates setting up and checking the original finite 
difference simul.it ion. It is, however, advisable to use double 
precision arithmetic in the coiiquiter code to minimire round off 
error . 

Soj^cMid, variable node spacing is permitted throughout. 
Although this feature is in the code published by dones and 
Pascoo (IP”!!, the simulation discussed below is generallv of 
a higher degree ap]M-ox imat i on . On the other hand, a simulation 
published by Ward et^ (19T’41, although of the same degiee 

approx imat i on as the one considered here, is tor a single 
valued constant nodal s]>acing. 

Thijd, the earth may be generally het erogeneous in two 
dimensions. The simulation is developed so that the nodes are 
placed at the corners of infinitely long two dimensional rec 
tangul.ar prisms. Each neighboring prism can be assigned a 
distinct value of resistivity. The gener.il finite difference 
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template then uses a five-point, second degree approximation 
for the most pathologic situation - the case of four dissimilar 
rectangular prisms abutting one another in the vicinity of the 
center point. Clearly, all cases of lesser pathology can be 
deduced from this situation and, in the case of a homogeneous 
material with equal nodal spacing, our general template reduces 
to the simple homogeneous finite difference simulation of the 
time-harmonic diffusion equation (Ward et al^. 1974). 

Fourth , and finally, machine coding using the general 
operator is enormously simplified because the boundary con- 
ditions to be applied between dissimilar materials are implic- 
itly handled by the relation. Therefore, to determine the 
coupling matrix, one needs only to specify the resistivity 
of each rectangular prism and then simply apply the operator 
sequentially centered at each nodal point. Although this latter 
feature is present in the code of Jones and Pascoe (1P"1), from 
their discussion one is not aware of the degree of accuracy of 
their approximation in the vicinity of discontinuities. Further 
more, even though Jones and Thompson (19'’4) acknowledge that 
the suggestions of Williamson e_t al^. (1 974) can improve the 
degree of approximation in a homogeneous region, they have not 
shown the degree of approximation in the neighborliood of dis- 
continuities. It seems appropriate therefore at the present 
state of literature to formulate a finite difference scheme with 
careful attention to the approximations involved. This is the 
intention of the following discussion. 



-25- 


OK 






Decoupling the Field Equations . It is v?ell known that 
Maxwell's equations in the M.K.S.A. system reduce to the fol- 
lowing forms: 


V X H = + iue)E (49) 

V X E * -iwy H (50) 

V • E « 0 (51) 

V • H = 0 (52) 


when one assumes the usual constitutive relations being 

the. true of d.c. conductivity), a time dependence of the form 
exp(iwt), and that e, the permittivity, and y, the permeability, 
are the same everywhere. Equation (49) leads to the condition 
that the total current is divergence-free, or that 

V • (J+ 3D/8t) = 0 (53) 

where J is the true conduction current and 3D/3t is the displace 
ment current. Assuming that the tangential msgnetic field is 
continuous across interfaces is sufficient to ensure that (53) 
is implicitly observed. A second boundary condition is that 
the tangential electric field is continuous across interfaces. 

We simplify the writing of (49) somewhat if we let o = ^ iwe 

where we call o the total conductivity and let it be equal to 
the sum of the d.c. conductivity, ^ , plus the term iwe. 

If one is in the principal coordinate system of a two-dimen 
sional inhomogeneity such that the structure strikes in the 


-26- 


OF POOR QUALITY 

y-direction, then for the class of source fields considered 
here, all derivatives of the fields with respect to y vanish. 
Moreover, the set of equations (49) through (52) d».;couples into 
two independent modes. 

One mode is determined by the set of equations relating 


and Hy through 


oE^ « -3Hy/3z 

(54a) 

oE = 3H /3x 
z y 

(54b) 

itoyH = -(3E /3z - 3E /3x) 
y X z 

(54c) 


and is knovvTi as the magnetic polarization, the H-parallel or the 
TM (transverse magnetic) mode. We see that if in some way the 
magnetic _■ ield component parallel to strike (H^) is determined 
everywhere, then and can be determined using (54a) and 
(54b) respectively. 

The other mode is determined by the set of equations re- 
lating H , H and E through 

X Z y 




- 

(3H^/3z - 

3H^/3x) 



(55a) 



iwpH^ 

= 3Ey/3z 




(55b) 



iw(jH^ 

= -3Ey/3x 




(55c) 

and 

is 

known as the 

electric 

polarization. 

the 

E-parallel 

or 

the 

TE 

(transverse 

electric) 

mode. We see. 

in 

this case, 

that 

if 

in 

some way the 

electric 

field component 

pai 

'allel to si 

trike 
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(E ) is determined everywhere, then H and H can be determined 
y X z 

using (55b) and (55c), respectively. 

Clearly the behavior of the fields for an arbitrary polariza- 
tion of the source field can be decomposed into two principal 
modes and, conversely, the two principal modes can be used to 
synthesize the solution for any arbitrary polarization. In the 
case considered here, the two-dimensional heterogeneity is 
oriented so that its strike is parallel to the direction of 
current flow in the source, hence we need to consider induction 
for the transverse electric (T5 or E-parallel) case. The TM 
(or transverse magnetic) case wil?. be considered in subsequent 
work. 

Operator for Electric Polarization . The electric polariza- 
tion mode is determined by the solution of equation set (55). 

We assume that the nodal point P is at the boundary between four 
electrically dissimilar regions as shown in Figure 6. In this 
case, we assume that in each of the regions the electric field, 

E^^, can be locally approximated by a quadratic function e^(x,z), 
where the subscript denotes the jth prism (j = 1 to 4) adjacent 
to the point P. Each quadratic function, 6j(x,z), can be ex- 
panded as a Taylor series in the vicinity of the point P, but 
in the jth region. For example, along either of the planes 
bounding region 1, 

e^(x,z) = e^ (0, 0) + x3ej^/8x + z3e^/3z 

(x^/2)3^ej/3x^ + (z^/2)3^e^/3z^ 


+ 


(56) 
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with similar expressions for e 2 » Cj and e^. Continuity of the 
tangential electric field requires that 


e^{0,0) = e^fO.O) = ejCO.O) = e^(O.O) = Cp , (57a) 

along with the following conditions at each of the nodes shown 
in Figure 7: 


e^(0,-dj^) 

- e^(0,-dj^) 

III 

(57b) 

1 

O 

= e2(-d^^,,0) 

- e^, , 

(57c) 

e.iO.dg) E 

ej(0,ds) = 

’ 

(57d) 

e^(dE,0) - 

e^(dp,0) £ 


(5'’e) 


Continuity of the tangent: al magnetic field across each interface 


implies, from (55b), that along the 1-2 interface 

ae,/9r = 9e^/9: . (58a) 

Along the 2-3 interface, (55c) implies 

9e.,/9x*9e^/9x . (58 b) 

Along the 3-4 interface, (55b) implies 

9e,/9z=9e^/9: . (5801 

And finally, along the 4-1 interface, (55c) implies 

9c, /9x - 9e,/9x . (5Si.n 

4 i 
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Expressions for each , written in the form (56), can be re- 
arranged to determine the appropriate first partial derivatives 
which are needed for substitution into relations (58a) through 
(58d). These Taylor series approximations can be combined and 
rearranged by making use of the fact that in each region j the 
second derivatives of the electric field are related through 
the diffusion equation 

3^e./3x^ + 3^e./3z^ = iuuo.e. . (59) 

3 3 3 3 

This leads to the following form for the electric polarization 
operator : 

(60) 

where 6g, the electric polarization induction parameter, is 
given by 

6^= io.y(Ojd^d^.* 0jd^4s-> cJjdsdg* a^d^dj)/4 . (61) 

Expression (60) is the finite difference operator we need for 
simulating the behavior of the electromagnetic fields at all 
nodal points interior to the boundaries of our modelling region, 
exclusive of the region of the source current itself. 
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Finite Difference Simulation of the Boundary Conditions, z « ± » 
It is straight-forward to replace the vertical derivatives 
in (46) or (48) with their appropriate finite forward or back- 
ward difference counterpart, such that they can be respectively 
written as 

Cp - (1/v.) (e^ - ep)/dg - 0 , (62) 

and 

However, these expressions suffer some loss of accuracy because 
the forward and backward difference formulae are accurate to 
second degree only at the mid-point of the layer which is 
adjacert to the nodal point P (Forsythe and Wasow, 1960). We 
can compensate for curvature effects on the field quantities 
through taking into account second derivative terms as shown 
in the following section. 



Top-most and Bottom-most Nodal Planes 


We assume that ey(x,z) can be represented by a Taylor's 
series expression in the medium immediately adjacent to and 
below a representa ".ive node, P, at the top-most nodal plane 
(z*z^) . The vertical derivative at can be expresseJ, to 

second degree accuracy, as 


3Ey/3z 


l(e^-ep)/d^l- (d,/2)3^ey/3z2 


(64) 


From the diffusion equation for E we have 


3^e^,/3z^ - iwua.e^ - 3^e^/3x^ 

> ^ y y 


(65) 


which can be substituted into (64) to yield, in terms of finite 
difference forms, 


(3E /3z) , = [(e -e )/d ] - 

> S p 5 


(ia)yo^d^/2)e^ + 




where the operator, 


+ (2/ [d^.(d^,+d^) ] }e^, , 

is the second finite difference simulation of 3‘‘Ey/3x" (Forsythe 
and Kasow, 1960). We now substitute (t>6) into (46) to obtain 
the following mixed boundary condition at z= 2 ^: 

Y^d^^e^ - [ (e^-e^) - iwyo^d^/2)e^ + (d^/2)D^^(ey)] •= 0 . (681 


An analogous condition for the bottom-most nodal plane can be 
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expressed as the following; 
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FINITE DIFFERENCE SIMULATION OF THE 
SOURCE FIELD WITHIN THE INTERIOR OF THE MESH 


We assume for the purpose of this discussion that the 
source is a thin horizontal current sheet such that, as des- 
cribed above, 


J 


y 


IpCosB 


X • 


(70) 


The geometry in Figure 8 represents a 5-point cell which con- 
tains the current sheet and whose center point, P, is immediately 
below the current sheet. We assume that for cells having dimen- 
sions small compared to the spatial wavelength of the source, 

Xg, one can reasonably approximate Ampere's Circuital Law in 

integral form 
* 

OH = 

with the expression 


J*dA 


(71) 


H ,d, - H ..d, 
xl 1 z2 2 


H -d- + 
x3 3 


H .d. 
z4 4 


I cosB d 
o XX 


(72) 


where H^^, for example with reference to Figure 8, represents 
the horizontal magnetic field component evaluated half-way along 
the line P-S and projected a.'.ong the length of the contour ele- 
ment dj^. Furthermore, we assume that the length d^ = (dj^,+d^)/2 
and that the respective magnetic field components evaluated 
half-way along each nodal spacing can be reasonably approximated 
by finite difference expressions in the form of centered dif- 
ferences. In other words, we express the left hand side of 
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(l/ia5y)(9Ey/37.) ^ (l/iwy)(3Ey/3x) ^(d5+dj^)/2 

-(l/iuy) (9Ey/9z) j^(d^+d£)/2- (l/iwy) (3Ey/3x) j^(d5+d^)/2 


« (1/iwy) [(d^+d£)/2] [(eg'ep)/^^] ♦ (l/iwy) [(d^+dj^)/2] (e£-ep)/dp] 
-(l/iu)y)[(d^+d£)/2]I(ep-ej^)/dj^] - (l/io)y) [ <:d 5 +djj)/ 2 ] [ (Cp-e^^.) . 

(73) 


The right hand side of expression (72) .-.an be integrated analyt- 
ically, 


cos [k (x-x^) ] dx , 


*p-V^ 


leading to 

J«dA -= (lp/k)[sin k (Xp-x^ + d^/ 2) - sin (Xp-x^-d^,/2) ] . 

(75) 

If k approaches 0, then the right hand side of (75) approaches 
(dui'^dc) /2 , as one would expect. 

Equating (73) to (75) provides the finite difference simu- 
lation of the source field condition at the nodal point P within 
the interior of the mesh. 
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Design Criteria for Mesh Spacing 

The case when the source field itself has a finite dimension 
imposes additional considerations for the problem of designing 
“dequate 2-D models. In general, the source field itself behaves 
according to the relation 

« exp (-kz) [cos (kx) + i sin (kx) ] , 

A 

where k» Zn/Xg, x is the horizontal distance and z is the ver- 
tical distance. Hence the mesh spacing, ^s, must be sufficiently 
dense to allow quadratic forms employed by the numerical algorithm 
to adequately represent trigenometric functions in the horizontal 
direction and exponentially damped functions in the vertical 
dimension. This suggests that to account for the finite dimen- 
sions of the source field, hs<\^/X2.S, for well-designed models. 

The induced field, of course, not only has a finii' wave 
number equivalent to that of the source but may also be mod- 
ulated by geometrical terms like 

exp (-d/6) 

where d is a distance parameter aud 6 is the depth of penetration. 
In some cases, close to electrical discontinuities the induced 
field may behave more like 1/r or 1/r . This often necessitates 
excessively close nodal spacing near electrical discontinuities. 

For the case of uniform source fields incident on a 2-D earth, 
one often increases the nodal spacing to very large values as 
one approaches the edges of the model. Gen. -illy, this is in- 
appropriate for finite source fields since, even a.‘ '-ne approaches 
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Testing the Alficrithm 

This section describes several of many tests to which the 
computer code was subjected. Figure 9 shows the behavior of 
only the source field in the interior of a rectilinear modelling 
region. The conductivity of the earth was set equal to the 
value of free*space and the source field was assigned a spatial 
wavelength of 200 km. The model values are shown as solid point 
in the figure and the straight line indicates the corresponding 
theoretical results. Besides showing the good agreement between 
the model and the theoretical results, this figure also shows 
the rapidity with which source fields having wavelengths of 
200 km or less are damped with distance below them - field 
amplitudes are reduced by an order of magnitude at a distance 
I'f onlv 7 0 km beneath the source. 

The effects of attenuation due to induction in an earth 
having a finite resistivity are shown in Figure 10. Mere tlu 
value of the electric field at a point within the earth is 
normalired by its amj'litude at the surface of a homogeneous 
half space having a resistivity of 1 ohm-m; the ]H'riod of the 
source is 1 sec and its siuitial wavelength is again 200 kn\. 

The values of tliese juirameters are chosen so that elect romagnet i 
induct ion effects clearl>’ dominate the effects of the finite 
spatial wavelength (i.e, [unio] k"l. The solid points cor 
respond to the model values and the straight line corresponds 
to tl'.e theoretical attenuation associated with a depth of 
penetration (61 of hlh' m. 



Indue t i on Near l.at or al He terogene ities 

The simple model shown in Figure 11 was used to examine 
the typical range of influence associated with anomalous fields 
near the vicinity of a lateral discontinuity excited by a finite 
source field. The source has a wavelength of SOO km, a period 
of l,e84 sec and a current strength of 1 amp/m. The resistivity 
on either side of the vertical contact is 50 ohm m and IS ohm m, 
respectively. Distances in Figure 11 are actually shov^^l to scale 
relative to the source wavelength and the relevant skin depths 
in each medium. 

Magnetic field values calculated for the model are shown as point 
in Figure IT; the open squares correspond to the values for 

the open circles correspond to the values for M, . The solid 
lines correspond to the theoretical and fields for the 
equivalent source field above a homogeneous half space having 
a resist ivitv of 50 ohm m; the vlashed lines correspond to ihe 
theoretical fields for a homogeneous half-s]’ace of IS ohm m. 

It should be clear that at some distance to the left of the 
contact, the numerical results for both horwontal and vertical 
magnetic field values fall on the cor respond ing theoretical 
curve for a 50 ohmm homogeneous half space. Similarly, at 
some distance to the right of tlie contact, the numerical re 
suits correspond to the results for a homogeneous half space 
of IS o hm m . 

What IS surprising is that the mod.el makes the transition 
between these two asvmptotic limits within a very small »listanct' 



of the contact. To all practical purposes, the transition oc- 
curs with 20 km or so of the contact. This implies that the 
scale length of the perturbation due to the contact is much 
less than either the corresponding depth of penetration in the 
representative media, or of the spatial wavelength of the source 
field. 

These results suggest that while fields may be perturbed 
in the close vicinity of lateral heterogeneities, at even modest 
offset distances from the contact the fields are more influenced 
by what is beneath the observer than by what is at some lateral 
distance away. This implies that in many cases, simple 1-D 
plane - layered models may be appropriately used in interpreting 
field data in a "piece-wise" sense, even when the total structure 
is more complicated laterally. 
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Magnetic Gradiometry Across Mega-Tectonic Feature s 

The fact that the fields induced by finite sources tend 
to decouple fairly rapidly with distance to either side of a 
vertical discontinuity prompted us to '•onsider the effects that 
this phenomenon has on magnetic gradiometric experiments. 
Schmucker (1970) and Kuckes (1973) have shown that one can 
determine an apparent conductivity (o ) for the earth through 
calculating the ratio of the horizontal gradient of the horizon- 
tal field (3H /3x) to the amplitude of the vertical field (H ) . 
The relevant expression when >> is 

(l/wy)l(3H^/3x)/H^|2 , 

We have applied this relation to calculating the apparent 
resistivity (p_ = l/o_) in the vicinity of a large scale discon- 
tinuity in the conductivity of the earth. Our model is simply 
a vertical discontinuity separating two quarter-spaces having 
resistivities of 30 ohm-m and 15 ohm-m respectively (Figure 13). 
A source field having a spatial wavelength of 20,000 km is 
located so that H is maximum over the contact and the gradient, 
3H^/3x, is essentially linear. The period of the source is 
10,800 sec (3 hr). After calculating the fields produced by 
such a model, two points were selected, point (1), in Figure 13, 
at -600 km, and point (2) at +400 km from the contact. A simple 
gradient, AH^/lO^m, was then calculated, providing an estimate 
for 3H^/3x. This r egional gradient was then used to calculate 
local response functions by employing first the vertica.1 field 


amplitude at point (1) and then the vertical field at point (2). 
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This allowed us to investigate the efficacy of employing local 
response functions in the vicinity of lateral discontinuities 
in the electrical conductivity. 

The results of our simple calculations are shown in Table 1, 
where the true resistivity beneath each site is compared with 
the magnetotelluric apparent resistivity and with the magnetic 
gradiometric apparent resistivity. The close agreement between 
these values is encouraging in that, at least in some cases, 
regional gradients of the horizontal field (even where the sites 
span a major discontinuity) can be combined with local values 
of the vertical field, to yield estimates of electromagnetic 
response functions. 
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Comparison of Satellite and Ground-Based Magnetic Variations 

Ke now consider the plane-earth equivalent of the situation 
illustrated in Figure 14. A satellite, at a nominal altitude 
of 400 km, passes above an ionospheric current system (at a 
height of 110 km), along a path which crosses a major geologic 
contact separating materials having a resistivity contrast 
of a factor of 2, not untypical of differences in gross con- 
tinental lithologies. We represent the surface of the earth 
by an infinite horizontal plane and consider only two-dimen- 
sional heterogeneities. Although our numerical model can handle 
more complicated geologic structures, here we consider only the 
case of a vertical contact separating two quarter- spaces having 
resistivities of 30 ohm-m and 15 ohm-m, respectively. 

Typical results from our model simulation for ionospheric 
currents flowing parallel to the geologic di scont inui ty are 
shown in Figure 15. Fields at the earth's surface are shown 
in Figure ISA, whereas fields at the satellite altitude are 
shown in Figure 15B. In both cases the results for a uniform 
source field are compared with those for a finite source field 

20,000 km) at a period of 10,800 s (3 hr.). 

At the earth's surface, there is little difference between 
the form of the H fields produced by the uniform and the finite 
source, respectively. The 1 fields produced bv cither source 
are also quite similar in form. Clearly the behavior of the 
fields near the vertical discontinuity are dominated by local 
induction - in other words, the depth of penetration dominates 


the scaling relation. 
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The major difference between the results produced by 
either type of source is that the field amplitudes associated 
with the finite source field are slightly biased to lower values 
due to the geometrical spreading factor exp(-kz). For a source 
with - 20,000 km, at a height of 110 km, the field amplitudes 
«t ground level will be reduced to 97 % of their maximum values. 
This is approximately the magnitude of the offset indicated 
in Figure 15. 

At satellite altitudes, the difference between the results 
for the uniform source field and those for the finite source 
field are more profound. The total horizontal field amplitude 
produced by the uniform source is approximately unity, suggesting 
thvit for remote sources in the magnetosphere the horizontal field 
is insensitive to the effects of lateral variations in the earth's 
conduct ivity . 

On the other hand, the total horizontal field amplitude 
produced by the finite source field is reduced by more than an 
order of magnitude. This is due to the fact that the satellite 
is above the source and, while the H-tield for the source itself 
is southward directed Ixl Pt satellite altitudes, the induced 
field is northward directed (,♦:<’). The two contributions tliercfore 
tend to cancel and the resultant H- field is quite small. At 
ground level the scale size of the anomalous H field is on the 
order of 300 km and api’'ears to be controlled by induction effects, 
whereas at satellite altitudes the scale size is on the order of 
IlPO 1200 Km and appears to be dominated bv the geometrical 
spreading factor cximK:!. 


-44- 


At satellite altitudes the vertical field produced by 
the uniform source field is quite small in amplitude and has 
a half-width of approximately 1000 km. At ground level the 
vertical field anomaly has a half-width of approximately 250 km 
and is approximately 5 times stronger in amplitude. 

The finite source produces an anomalous vertical field at 
satellite altitudes which is of the same order of magnitude as 
that produced by the uniform source, the difference being an 
apparent skewing of the sTiape of the finite 2 field due to the 
fact that the anomalous field locally associated with the ver- 
tical discontinuity is superimposed on a broad regional field 
which varies sinusoidally. This regional field is positive 
for x<0 and negative for x>0. The local anomalous field is 
negative, hence the sum appears to be skewed to the right (i.e. 
to X > 0) rather than having the symmetrical form of the perturba- 
tion associated with the uniform source. 
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TABLE 1. Comparison of Local Response 
Parameters Using Regional 
Horizontal Field Gradients 

Source Field Wave-Length ■ 20,000 km 

Period ■ 10,800 sec (3 hr) 


Resistivity, ohm-m 

Magnetic 

True Magnetote 1 lu . ic Gradiometr ic 


Site 1 

30 

29 

27 

Site 2 

15 

15 

16 
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figure CAPTIONS 

Figure 1. Geometry of the problem considered in this paper. 

An equivalent current source, Jy(x) flows at some height 
parallel to the earth's surface and parallel to the axis 
of a 2-D inhomogeneity which has a characteristic scale 
dimension L. 

Figure 2. One approach to the problem when typical wavelengths 
for the source field, X^, are much greater than the dimen- 
sions of the inhomogeneity L. Values for the source 
field and some average reflectance for the earth are 
assumed along the top plane of the modelling region. 

Figure 3. A situation where only the source field is assumed 
to have appreciable magnitudes in some region of space 
(or the modelling region) and where all fields are assume; 
to tend to iero at the remote edges of the model. In 
general this requires that << L. 

Figui 4. Details of the model considered in present paper. 

ht finite source is located within the modelling region 
and characteristic impedance conditions for a finite 
source are imposed at the top-most and bottom-most edges. 
Earth structures, while having to be plane- layered , may 
be different at the left and right edges. 

Figure 5. An equivalent current shec*^ source flowing parallel 
to the earth's surface. 
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Q> ^4. - 

Figure t>. The situation of four dissimilar materials abutting 
at a contact. 

Figure 7. Convention for indicating nearest neighbor nodes 
surrounding the central node P (following Forsythe and 
Wasow; 1960). The various media surrounding the central 
node P are indicated by j = 1 to 4. 

Figure 8. Geometry and conventions used to simulate a dis- 
continuity in the horizontal sheet current densitv J (xl , 
flowing above, but adjacent to a nodal point P. The 
shaded region represents the amplitude of an infinitesimal 
thin sheet current as a function of z. The dashed 
contour with arrows represents the line integral used to 
evaluate Ampere's Law. 

Figure 9. Comparison of model results (solid points') and 

analytical results (solid line) for the source field alone. 

Figure 10. Comparison of model results (solid points) and 

analytical results (solid line) when attenuation of the 
field due to electromagnetic induction dominates its 
behav i or . 

Figure 11. Model geometry used to illustrate the coupling of 
fields across a contact. 

Figure 12. The beha\ior of model field va.ues (sec Figure 11) 

for (seynares) and il„ (circles) across a contact cominired 
with analyticai results for a homogeneous earth having a 
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resistivity of 30 ohm-m (solid line) and 15 ohm-m (dashed 
line) . 

Figure 13. A simulation of magnetic gradiometry across a ver- 
tical contact separating 30 ohm-m and 15 ohm-m materials. 
The source field has a spatial wavelength of 20,000 km 
and a period of 10,800 sec (3 hr). Results are summarized 
in Table 1. 

Figure 14. Geometry representing a satellite orbiting at a 
height of 400 km, an ionospheric current source at a 
height of 110 km, and a vertical electrical contact 
separating gross earth lithologies having resistivities 
of 30 ohm-m and 15 ohm-m, respectively. 

Figure 15. Comparison of simulated satellite and ground-based 
observations of magnetic variations. The H field is 
horizontal and directed toward increasing x; the Z field 
is vertical and directed downward. The signs in brackets 
are associated with the various branches of the curves 
used to denote the amplitude of the respective field 
amplitudes. Panel A: H and Z fields observed on the 

earth's surface; Panel B; H and Z fields observed at 
satellite altitudes (400 km). All field amplitudes rep- 
resent the sum of the source and the induced field com- 


ponents. 
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TABLE 1. Comparison of Local Response 
Parameters Using Regional 
Horizontal Field Gradients 

Source Field Wave-Length * 20,000 km 

Period » 10,800 sec (3 hr) 

Resis* ^ '"ity , ohm-m 




Magnet ic 

True 

Magnetotelluric 

Gradiometr ic 

Site 1 30 

29 

27 

Site 2 15 

15 

16 
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